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We consider the cumulant expansion of the PAM employing the hybridization as perturbation 
(Phys. Rev. B 50, 17933 (1994)), and we obtain formally exact one-electron Green's functions 
(GF). These GF contain effective cumulants that are as difficult to calculate as the original GF, 
^ ■ and the Atomic Approximation consists in substituting the effective cumulants by the ones that 

' correspond to the atomic case, namely by taking a conduction band of zeroth width and local 

hybridization. This approximation has already been used for the case of infinite electronic repulsion 
U (Phys. Rev. B 62, 7882 (2000)), and here we extend the treatment to the case of finite U. The 
method can also be applied to the single impurity Anderson model (SI AM), and we give explicit 
expressions of the approximate GF both for the PAM and the SIAM. 
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I. INTRODUCTION 



In this work we discuss approximate Green's Functions (GF) for the Periodic Anderson Model (PAM), obtained 
by starting from a formally exact expression and approximating a component of this expression by the corresponding 
^ exact solution of the atomic problem. We have already employed this technique in the limit of infinite repulsion U of 

• the localized electrons and here we shall extend the technique to the case of a finite U. We call this technique 

the Atomic Approximation, not to be confused with the atomic solution of the problem. 
The Hamiltonian for the PAM is 



k.o- j,a 
+ ^ X] + Hh, (1) 



> 

where the operators ^ and Ck,cr are the creation and destruction operators of conduction band electrons (c-electrons) 

\ with wave vector k, component of spin a and energies i^k.u- The /J^ and fj^„ are the corresponding operators for 
• ■ the /-electrons in the Wannier localized state at site j , with spin component a and site independent energy E^. The 
, third term is the Coulomb repulsion between the localized electrons at each site where rij^a = fjcrfj,^ is the number 



of /-electrons with spin component a at site j and the symbol a denotes the spin component opposite to a. The 



' fourth term Hh describes the hybridization between the localized and conduction electrons 

. si : Hh=J2 iVjM,<rflC,.,<r + V*^^,Clj,,.), (2) 

H ' . 

" with a coupling strength given by 



VjM,a = ^^K(k) exp (ik.R,), (3) 

where A^^ is the number of sites in the system and V'o-(k) is independent of the wave vector k when the mixing is 
purely local. 

If we consider that the local repulsion between /-electrons is infinite {U — > oo), so that the double occupancy at any 
site is zero, we can employ Hubbard X operators to make disappear the term proportional to U in the Hamiltonian. 
To this purpose we consider first the definition of the X operators: the Xj^ha transforms the state | a > at site j into 
the state | 6 > at the same site, and we assume that | a > and | b > are eigenstates of the number of electrons. We say 
that Xj^ba is of the Fermi type when | a > and | b > differ by an odd number of Fermions, and that it is of the Bose 
type when they differ by an even number of Fermions. By definition, two AT-operators of the Fermi type at different 
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sites anti-commute, and commute when at least one of them is of the Bose type. The algebra of these operators when 
they are at the same site is defined by their product rule 

^j,ab ■ Xj,cd — Sb,c ■ Xj,ad, (4) 

and they are neither Fermions nor Bosons. For infinite J7, the only /-electron states at any site j are the vacuum 
I j, > and the two states | j, cr > that have one electron with spin component cr, and the only Fermi type operators 
that we shall need in this case are Xj ^a and their Hermitian conjugates Xjuo = -^j oa- Projecting H into the subspace 
without doubly occupied f-electron states we obtain the PAM Hamiltonian for infinite U: 

H — ^ -Ek.trC'^^g.Clccr + ^ Ej^Xj^^a- 

+ E (^.,k,.4o.t^k. + ^tk.^cL^xo.) (5) 



where Xj-.o-o- = oa-^j,o<j is the projector into the state \ j,a > . The identity relation in the reduced space of the 
localized states at site j is 

Xj,oo + Xj^ca + = I (6) 

and its statistical average gives the conservation of probability in that space of states. 

The generalization of Eq. (O to the case of several configurations with a rather arbitrary choice of states is 

H = ^ E]s_,aC\^^^C\^,a + ^ EjaXj,aa + Hh = Ho+Hh, (7) 
k,(T 

where ""^ 

jba,k(T 

The a and b summations are over all the states | a > and | b > that we want to include in the model, and the only 
restriction is that any hybridization constant must vanish unless state | a > has just one electron more than the state 
I 6 >: this last condition is necessary to satisfy the conservation of electrons. In this general case, the energies Ej^a 
include all the Coulomb repulsions of the type described by the third term in Eq. ([T|). 

To abbreviate, we can write Eq. ([5]) in the interaction picture in the more compact form: 

H^ir)^J2^{l,l')Yil)Y{l'). (9) 

1,1' 

where 

Yil) = y^(T) = cxp {tHo)Y^ exp (-ri/„) (10) 

is the operator Y^ in the interaction picture (the subindex 7 is discussed in more detail after Eq. (jl5p ). The only non- 
zero coupling coefficients Vil, I') are those that correspond to the correct combination of indices I and I' in Eq. ^ and 
a factor 1/2 is not necessary in Eq. ([5]) if we choose to retain only terms in which Y{1) corresponds to the /-electrons 



^ For Eq. JTJ with finite U there are four states |0}, |+), |— ) ,\d) = \-\ — ), and we have fa = Xq„- + crX^jj. { cr = 1 corresponds to + and 
(T = —1 to — ). When we write Eq. [2] employing Xoo- and X^d it appears aX^d rather than X^d, and we have to write the Vj,ba,]i,tT in 
Eq. [8] in the following way: V,,<^d,k,<T = o-V'j,OCT,k,<T = o"V,,k,CT- 
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and Y{1') to the conduction electron (to achieve this ordering in the second term of the parenthesis in Eq. ([8]) one must 
anti commute two Fermi type operators, and the corresponding minus sign is absorbed into a redefined hybridization 
constant). 

As we are interested in the Grand Canonical Ensemble of electrons, we should replace the total Hamiltonian H by 



n = H-^JYl ^l.^K- + ''-^i-- >^^o + Hh (11) 

[ k,CT ja J 

where Xj aa is the occupation number operator of state | a > at site j, and Va is the number of electrons in that state. 
This transformation is easily performed by changing the energies Ej^a of all ionic states | a > into 

^j,a = Ej^a - M^^a (12) 

and the energies i?k,cr of the conduction electrons into 

e{-k,a)^Eu,a-l^ (13) 

II. THE GREEN'S FUNCTIONS IN IMAGINARY FREQUENCY FOR SEVERAL PARTICLES. 

In this work we shall consider Green's Functions ((GF) of imaginary times of both conduction electron operators 
Cko- and Hubbard operators Xj ^a, and a general GF can be written employing the operators: 

^(7i,ri;--- ;7n,r„) = <^(f(7i,ri)...y(7„,T„))^^ , (14) 

where 

y (7, t) = exp {tH)Y^ exp {-tU) (15) 

is defined for /3 > r > 0. Besides the Fermi- like operators Ky that appear in Hh, we shall also consider Bose-like 
Hubbard operators that do not change the number of electrons. At this point it is necessary to be more specific about 
the argument 7 of the operators Ky in Eqs. (|10p and ([15]). When the corresponding Y^ is a Fermi type Xj.ha, we use 
7 = (/; j, Oi, u), with u — —, and the single index a identifies the transition \ a > \ b >, with the same restriction 
stated after Eq. ([5]), namely that state | a > has just one electron more than the state \ b >. The inverse transition 
(operator A'j^^) is described by the same a but with u = +. The j identifies the site, r is the imaginary time (cf. 
Eq. pUj) and Eq. (HH])) , and / is only used when necessary to avoid confusion. When Yj is Cka we use 7 — (c; k, a, u) 
with u = — and change to m = + for C^^.. It is not necessary to assign a u parameter to the Bose-type operators, but 
to unify the notation we shall keep the u and put always u = 1 for these operators. The only restriction on the two 
states I a > and | 6 > of the transition a = (6, a) for Bose type operators, is that they should have the same number 
of electrons. 

One can not use Feynman type expansions for the GF in Eq. (jl4p . because the Hubbard operators are not Fermi 
operators, and we shall use a cumulant expansion ^3^ that is an extension of the one derived by Hubbard for 
his model. The diagrammatic expansion of the GF is obtained employing the Theorem 3.3 from Reference @, that 
expresses the GF as the sum of the contributions of all the topologically distinct and vacuum free graphs, drawn 
according to Rule 3.4 of that reference. The corresponding contributions are calculated with Rule 3.6 of [s!], and in 
this section we shall summarize some details of these GF calculation. 

To avoid repeating the same term in Eq. (|9]) we assumed that V{1, V) is non-zero only when the first index corresponds 
to an A"-operator. These coefficients do not depend on r or t', and to abbreviate it is convenient to introduce 
a, k, CT, u) in Eq. ^ : 

v{j,a,'k,a,+) = V{f;j,a,+;c;'k,(T,~) = Vj,ba,Ka, 

v{j,a,k,a,~) = -VU;j,a,-,c;k,a,+)^V*- - ^'^^ 



,6a,k,(7 ■ 
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The minus sign that should multiply into V*^ - , because we anti-commuted two Fermi-type operators from Eq. ([8]) 

in the corresponding terms of Eq. will be absorbed in the rules for the sign of the graph contributions when ^ = 
(cf. Appendix El 2. 

To Fourier transform with respect to time the GF of Eq. p^ . it is essential that they obey the boundary condition 

^ (y (71, Ti) • • • y(7„ r, = /3) • • • y(7„, ^n)) = 

± ^ (f(7i, Ti) • • • y(7„ r, = 0) • • • f(7n, ^»)) (17) 

with respect to all the operators i^(7i, ti) • • • Y{jn, Tn), where the —{+) corresponds to Fermi-like (Bose-like) operators 

When Eq. (|17p is satisfied for all the variables and H does not depend on r, we can treat the GF as periodic (anti 
periodic) with period /3 in r, for all Bose-like (Fermi-like) operators K(7,r), and we then write 



(y(7i,Ti)...f(7„,T„)) 



X exp [-i(wiTi H h w„T„)] (18) 

The frequencies ujj are different for the two type of operators Ky: 

TTVj , f jy 0, =f2, t4 • ■ ■ Bose-like 

Lo.j = — where < i o c • n (19) 

^ f3 [ = 1, =F3, t5 • ■ ■ Fermi-like. ^ ^ 

The notation of the Fourier coefficients in Eq. ()18p is purely symbolic, because the r-ordering (. . .)+ has no meaning 
there. 

A. Rules for reciprocal space and imaginary frequencies 

To Fourier transform the spatial dependence one has to remember that the c-operators are already in reciprocal 
space, so it is only necessary to transform the /-operators. For a GF with r operators of the /-type (Fermi-like or 
Bose-like) and n — r operators of the c-type we write in an abbreviated notation 

(y(/, r; 1) • • • y(/, r; r)Y{c, r; r + 1) • • • f(c, r; n)) - 
P^^Ns^ E E exp[— i(kiWiRi + • • • + kpU^Rr) - ^(wiTi -I- • • • -I- a;„T„)] 

X ^ (y (/, Lo;l)--- Y{f, uj; r)Y{c, u-r + !)■■■ F (c, uj;n)) ^ ^, (20) 



^ Note also that a factor (—1)" appears in the perturbation expansion contribution of any graph of order n, i.e. with n internal edges 
(cf. the cumulant expansion for the Ising model in Ref. 5], where this sign has been included in the interaction constant in its Eq. 
(2)) We have then added a factor (—1) to every internal edge, and therefore this extra factors would only change the sign of a graph's 
contribution when it is of odd order This sign appears explicitly in the expansion of the PAM in Q (cf. Eqs. (3. 8), (3. 11) of that 
reference) but it was left out from the diagrams contribution by an oversight. Note that this sign does not depend on the Fcrmionic 
character of the X operators. 
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where Rg is the position of site js, Y{f,T;s) = Y{f;js,as,Us,Ts), Y{c,t;s) = y(c; k^, cr^, Us, Ts), and we substitute 
the Ts by Us in Y(f, w; s) and Y{c, uj; s), as well as jg by in Y{f, uj; s). With the same notation, the inverse relation 
is then 



^ uj;l)--- Y{f, UJ- r)y(c, c^; r + 1) • • • y(c, n)) ^ = 

-r 

j3^^ Ns ^ / dri - ■ ■ dTn exp[+i(kiwi • Ri + • • • + k,.Mr • Rr) + j(wiTi + • • • + iuj„t„)] 

X ^ r; 1) • • • r; r)f(c, r; r + 1) • • • y(c, r; n)) . (21) 

The present definition is slightly different from Hubbard's because we include the parameter u = ±1 into the 
spatial part of the exponential in Eqs. (PO)) and (PTt . 

From the invariance under time translation (i.e. H does not depend on r) one can show that the GF in Eq. (PT|) 
vanishes unless 



wi+waH hw„ = 0. (22) 

To prove the corresponding property for the wave vectors kj in Eq. (|2ip. it is necessary to transform first the 
c-operators into the Wannier representation 



^exp(-zk.R,)CL. 



(23) 



Substituting Eq. (|23|) (or its Hermitian conjugate for Cjcr ) into the GF in the r.h.s. of Eq. (|2T|) . and employing the 
invariance under lattice translation, one finds that the GF in Eq. ([?!]) vanishes unless 



kiui + k2U2 H h k„it„ = 0. (24) 

It is clear that the relations in Eqs. (|20ll2ip can also be employed for the corresponding cumulant averages. When 
Hh = many Y^ are statistically independent, and the only cumulants left in the imaginary-time and real-space 
expansion (cf. Rule 3.6 in reference [J]) must either have all their Y^ of the /-type and at the same site, or else 
have all of the c-type with the same k (and same a when Ho is spin independent). Because of the invariance of 
the system under lattice translations, the local cumulants that appear in Rule 3.6a in f§\ are independent of the site 
position and it is not necessary to take their spatial Fourier transform; on the other hand, the Y^ of the C-electron 
cumulants of Rule 3.6 b' in Q have been already transformed. From these two facts it follows that to obtain the 
Fourier transformed version (in reciprocal space and imaginary time) of Rule 3.6 in [sj it would be sufficient to apply 
only the transformation from time to frequency (cf. Eq. (1181) ) to the cumulants in that rule. 
To set the notation we write 

{Y{f;j: ap,Up,Tp) ■ ■ ■Y{f;j,ai,ui,Ti))^^^ = 
/3"5 ^ exp[-i{ujiTi -\ \-uJpTp)] 

{Y{i, ap,Up,ujp)--- Y{j, ai , ui, wi))_^^^ (25) 

and 



^ The parameter u was defined after Eq. I I15I I. as well as in the sentences after Rule 3.4 in reference Q. The u was convenient to organize 
our calculation, but we did not use it in the temporal part of the exponential because it was not particularly useful there. 
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((r(c; k, (J2,~U2, T2)Y{c; k, ai,ui,Ti))^) ^ 
^ exp[-i(wiTi +UJ2T2)] 



((C(k, a2,-U2,uj2)C{k, cri, mi, wi))^)^ 



(26) 



Note that the invariance under time translation guarantees that Eq. (|22l) would be satisfied for the frequency 
dependent cumulants of Eqs. (p5)) and (p6)) . To proceed with the transformation of Rule 3.6 in [3], we use the 
prescriptions summarized above to express the GF in the r.h.s. of Eq. ()21|) as a sum of terms, each corresponding to 
the contribution of some graph. In each term one introduces Eqs. (|25p and (|26|) and then performs explicitly all the 
integrations over r and all the non- restricted summations over the sites j. 

In each integration over t there are two possibilities: the r corresponds either to an external operator or else to an 
internal line. When the Tj corresponds to an external operator Y{'yj,Tj), the Eq. pT|) provides the integration, and 
the integrand has two factors: one eicp{iujjTj) from Eq. (j2ip and another exp(— iw^Tj) from applying Eqs. (|25p and 
(|26p to the cumulant of Rule 3.6 that contains the external operator Y^jjjTj). As both ujj and Ws are of the same 
type (cf. Eq. (fT9)) ). the integral vanishes unless cjj = w^, and from the sum over all the cug in Eqs. ((25)l and (|26|) only 
the external frequency cdj remains. 

When the belongs to an internal line, the integration comes from the perturbation expansion (cf. Eq. (IBip ). 
and the integrand is exjp[—i(ujs + Ws)ts] where uJs and ujg come from expanding with Eq. (j25l) or Eq. (j26l) the two 
cumulants of Rule 3.6 in 0] that contain the C-operator and the X-operator of the internal line. The integration is 
again zero unless ujs + = 0, and one can then associate only one of these two frequencies to the internal line in the 
transformed rules. 

In Eq. pT|) we have applied the spatial Fourier transformation to the external X-operators, which together with 
items (3a) and (2d) of Rules 3.5-3.6 from Q imply a sum over all the sites in the lattice. It is then convenient to 
write explicitly the dependence with Rj of the coupling constants of Eq. ([TBI) • 



and one then obtains the following Rule. 
Rule 3.7 

To calculate the contribution of any diagram obtained from Rule 3.4 of Q 

1. Assign to each internal line a momentum kg, a frequency ujs, a spin as and an index Us- Assign dummy labels 

as and ±Us to the X-operators at the FV side of the internal line, and dummy labels k^, as and ^Us to the 
C-operators at the CV side. Use +Us and +ujs at the side of the edge to which points the arrow (cf. item iv of 
Rule 3.4 in Q) and —Us and —Ug to the opposite side. 

Assign to the external lines the labels of the corresponding external operators, namely the momentum k^, 
frequency lOs, index Ug and also the transition as = (fos,as) for X-operators and the spin component as for the 
C-operators (we use always +Ur and +uJr for the external lines). 

2. Form the product of the following factors: 

(a) For each FV with lines s = 1, 2, • • • ,p running to that vertex (both internal and external) the factor'' 



where ks, Ws, as and Us are the momentum, frequency transition, as = (6s, a^) and index Us labels of the 
X-operators associated to line s (always +Us and +uJs for the external lines). 

(b) For each CV a factor 



v{j, a, k, a, u) — V{a, k, a, u)Ns ^ exp(iuk.R_,), 



(27) 



NsA{±Upkp ± • • • ± U2k2 ± uiki) 




(28) 



To simplify the notation wc use A (x) = when x but A (0) = 1. 
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((C(k'i,a'i,-<,-c^'i) C(ki,ai,ui,c^i))+)^, (29) 

where ki, cti, ui and uji are the parameters of the edge with the arrow pomting towards the CV. As we 
discussed before, this cumulant vanishes unless ki = k'j^, ui = u\ and wi = w'^. When the Bloch states | k, > 
are eigenstates of "Ho, we have also ai — a\ and the factor above (cf. footnote |8] in Appendix |F]) is equal to 



1 

zwi + uie(ki, cTi) 

X <5(ki,k'i)5(ui,<)<5(ai,a'i)(5K,c^'i) (30) 

where the parameters with sub index 1 correspond again to the edge with the arrow pointing towards the CV 
(when the outgoing line is external with given u and w, we put — u'j^ = u and — a;). 

(c) A factor (—1) F(a, fc, ct, ±ti) for each internal lin(Pwith labels a, ±m at the FV site and labels fc, u and =Fu at 
the CV side, as written in f3| (cf. Eq. . 

(d) A factor ±1 determined by the rules in Appendix C. 

(e) A factor l/g determined by the rules in Appendix D. 

(f) A factor Xj^fWa for each external line running to a FV. 
3. Sum the resulting product with respect to 

(a) The momenta k^, the frequencies and the indices Us of all the internal edges. Divide each sum over momenta 
into -JWs. 

(b) The labels as of the AT-operators at the FV side of all internal lines. 

(c) The label of the C-operators at the CV side of all internal lines.* 

Two points should be stressed: i) The frequencies of each local cumulant in 2. a satisfy Eq. p2|). thus reducing by 
one the number of frequency summations at each FV. ii) The rules are also valid for vacuum graphs, and are employed 
to calculate the GPF with the Linked Cluster Theorem. 



B. Rules for real space and imaginary frequencies (Valid for impurities) 

We shall transform Fourier the imaginary times of the diagrammatic expansion calculated with Rule 3.6 in but 
leave the real space description of the local sites unchanged. 

We employ the Rule 3.4 in fs*] for drawing the nth-order graphs for the cumulant expansion. The following relations 
give the Fourier transforms, following the same definitions employed in pi] 



^ (y(/, r; 1) • • • f(/, r; r)f (c, r; r + 1) • • • Y(c, r; ")) 
^ exp[-i(a;iri H |-WnT„)] 



UJ\ ■■■UJn 



X ^ (y (/, 0.; 1) . . . y (/, cj; r)y (c, Lu;r + 1)--- Y{c, c.; n)) , (31) 

where Y{f,T;s) = Y{f;js,as,Us,Ts) , Y{c,t;s) = Y{c;ks,as,Us,Ts), and we substitute the r by w in Y{f,uj;s) and 
Y{c,uj;s), but keep the js in Y{f,u!;s) because here we do not transform the GF into reciprocal space. With the 
same notation, the inverse relation is then 



^ Note that in Reference 01 we use V{a,k,(T,±u) in reciprocal space (rule 3.7, item 2.c) and v{a, k, a, ±u) in real space (rule 3.5, item 
2.c) 
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^ uj;l)--- Y{f, u; r)Y{c, w; r + 1) • • • Y{c, u; n)) 

3'^ I dri • • • / dT„ exp[i(wiTi H h w„t„)] 

Jo Jo 



X ^(f(/,T;l)...f(/,r;r)y(c,T;r + l)...y(/,r;n))^^ . (32) 

From the invariance under time translation (i.e. % does not depend on t) one can show again that the GF in 
Eq. ([5^ vanishes unless Eq. ([^ is satisfied (i.e. cji + a;2 + • • • + a;„ = 0). 

It is clear that the relations in Eqs. PHI32I) can also be employed for the corresponding cumulant averages. When 
Hh — many are statistically independent, and the only cumulants left in Rule 3.6 must either have all their Y^ of 
the /-type and at the same site, or else have all their Y^ of the c-type with the same k (and same a when T-Lo is spin 
independent). We shall then apply the transformation from time to frequency (cf. Eq. (3.29) in ^) to the cumulants 
in that rule. To set the notation we write 

(^{Y{f;j, ap,Up,Tp) ■ ■ ■Y{f;j,ai,ui,Ti))^ ^ = 
/3"5 ^ exp[-i(iLJiri H V tOpTp)] 

(^{X{j,ap,Up,ujp) ■ • •X(j,Qi,ui,cji))^^^ (33) 

and 

(^(y{c; k,a2,--U2,T2)Y{c; fc, cti, mi, n)^ = 

l3~^ ^ exp[~i{u>iTi +LJ2T2)] 
LJ1UJ2 

C{k,a2,~U2,u}2)C{k,ai,ui,uji)j \ (34) 



Note that the invariance under time translation guarantees that Eq. (|22l) would be satisfied for the frequency 
dependent cumulants of Eqs. (j33p and p4p. To proceed with the transformation of Rule 3.6, we use the prescription 
given in Rule 3.4 to express the GF in the r.h.s. of Eq. (|32l) as a sum of terms, each corresponding to some graph. 
In the contribution of each graph one introduces Eqs. ((33)) and (p4)) and then performs explicitly all the integrations 
over r while the non-restricted summations over the sites j remain expressed formally. In each integration over r 
there are two possibilities: the r corresponds either to an external operator or else to an internal line. When the 
Tj corresponds to an external operator Y^-jjjTj), the Eq. ((32)) provides the integration, and the integrand has two 
factors: one exp{iu}jTj) from Eq. ((5^ and another exp(— iw^Tj) from applying Eqs. ((55)) and ((Ml) to the cumulant of 
Rule 3.6 that contains the external operator Y{'jj,Tj). As both LUj and Wg are of the same type (cf. Eq. (3.30) in 
Q), the integral vanishes unless ujj = ujg, and from the sum over all the oJs in Eqs. ()33|) and (1341) only the external 
frequency ujj remains. When the Ts belongs to an internal line, the integration comes from the perturbation expansion 
(cf. Eq. (3.11) in Q), and the integrand is exp[—i{ujs -|- a;g)r5] where ujs and u)'^ come from expanding with Eq. p3p 
or Eq. ((Mj) the two cumulants of Rule 3.6 that contain the C-operator and the X-operator of the internal line. The 
integration is again zero unless cj^ + cjI, = 0, and one can then associate only one of these two frequencies to the 
internal line in the transformed rules. Before explicitly giving the rules for the GF calculation, it is convenient to 
recall the definition Eq. ((TC)) of the coefficients v{j, a, k, a, u): 

v{j,a,k,(T,+) = V{f;j,a,+;c;k,a,-) = V^^^^i:^^, 
v{j,a,k,(T,-) = -~V{f;j,a,-;c;k,a,+) = V*^^^^^^. 

The minus sign that should appear with V*^ - because we anti-commutcd two Fermi-type operators from Eq. (2.8) 

in 0, will be absorbed in the rules for the sign of the graph contributions when ^ = (cf. Appendix)^. We can now 
give the rules without further discussion. 

Rule 3.7a To calculate the contribution of a diagram obtained from Rule 3.4 
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1. Assign to each FV a site label js- To each internal line a momentum kg, a frequency utg and an index ztug- 

Assign dummy labels and ztu^ to the X-operators at the FV side of the internal line, and dummy labels k^,, 
(Ts and =FUs to the C-operators at the CV side. Use +Us and +ujs at the side of the edge to which points the 
arrow (cf. item iv of Rule 3.4) and —Ug and —LUg to the opposite side. Assign to the external lines the labels 
of the corresponding external operators, namely the frequency LUg, index Ug and either the site jg and transition 
ag — {bg,ag) for X-operators or the momentum and the spin component ag for the C-operators (we use 
always +Ur and +ujr for the external lines). 

2. Form the product of the following factors: 

(a) For each FV with lines s — 1,2, ■ ■ ■ ,p running to that vertex (both internal and external) the factor 

{{X{jp, ap, ±Up, ±ujp) ■ ■ ■ X{ji,ai,±ui, , (35) 

where jg, ujg, ag and Ug are the site, frequency transition, ag = (bg,ag) and index Ug labels of the X-operators 
associated to line s (always +Ug and +LUg for the external lines). 

(b) For each CV a factor 

{iCiK,a[,~u[,-Lu[)C{k,,auUi,LUi))^)^, (36) 

where ki, cti, ui and uji are the parameters of the edge with the arrow pointing towards the CV. As we 
discussed before, this cumulant vanishes unless ki = k'j^, ui = u'l and uii = lu[. When the Bloch states | k, cr > 
are eigenstates of Ho, we have also CTi = (t[ and the factor above (cf. footnote |S] in Appendix |F]) is equal to 

1 

iuji + uie(ki, CTi) 

X 5{ki,k[)S{ui,u[)5{<7i,a[)S{uji,uj[) (37) 

where the parameters with sub index 1 correspond again to the edge with the arrow pointing towards the CV 
(when the outgoing line is external with given u and uj, we put —Ui = u and — uj). 

(c) A factor (—1) v{j,a,'k,a,ztu) for each internal edg^ joining a FV at site j with labels a, ±u to a CV with 
labels k, a and ^u. 

(d) A S{js,ji) for each external line X-operator at site ji running to an FV site labeled with js. The labels jg are 
dummy labels, but the Kroenecker deltas in the present item take care of fixing its value when there is an external 
line running to a FV. 

(e) A factor ±1 determined by the rules in Appendix C. 

(f) A factor 1/g determined by the rules in Appendix D. 

3. Sum the resulting product with respect to 

(a) The site labels jg of all the FV. 

(b) The momenta k^, the frequencies ujg and the indices Ug of all the internal edges. 

(c) The labels ag of the X-operators at the FV side of all internal lines. 

(d) The label ag of the C-operators at the CV side of all internal lines.* 

Two points should be stressed: i) The frequencies of each local cumulant in 2. a satisfy Eq. ([22l) . thus reducing by 
one the number of frequency summations at each FV. ii) The rules are also valid for vacuum graphs, and are employed 
to calculate the GCP with the Linked Cluster Theorem. 
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III. THE EFFECTIVE CUMULANT 



The general GF in reciprocal space and imaginary frequency that we shall use is 
^/■''■^(k, a, u, w; k', a', w', w') = 

^ E e'^P ("k • R + uV ■ R')] dT J^^ dr' exp [i {ujt + Jr')] (9 (/; j, a, u, r) Y (/; /, a', w', r')) 



(38) 

where R' is the position of the site /, and in particular we need the transform of (^(^Xj,a{Ti)X^^, ^,{0)) ^ , i.e.: 



Q^-f {k^a^u = — , cj; k', a', w' — +,a;') . Employing Eqs. ( I22I24I) and the conservation of the number of electrons, we 
can abbreviate 

^■^■^(k, a,u = -,ij; k', a', u\ w') = (k, iuj,) A + u') A (uk + w'k') A {uj, + ^) (39) 

where uj and uj' are given by Eq. (|19p . 

In the calculation with the usual Fermi or Bose operators, the one-particle propagator of the f-elcctron is given by 
a sum of diagrams of the type shown^ in figure [1^ but with each local vertex corresponding to the sum of all "proper" 
(or irreducible) diagrams [sQ- The same result is found in the cumulant expansion of the Hubbard model for d — > oo 

m when the electron hopping is employed as perturbation. The vertices then represent an "effective cumulant" 
M2^J {zn), that is independent of k because only diagrams of a special type contribute to this quantity for d — >■ oo. 

In the cumulant expansion of the Anderson lattice [3|] we employ the hybridization rather than the hopping as a 
perturbation, and the exact solution of the conduction electrons problem in the absence of hybridization is part of 
the zeroth order Hamiltonian. For this reason it became necessary to extend Metzner's derivation to the Anderson 
lattice for U — oo, and we have shown[10] that the same type of results obtained by Metzner are also valid for this 
model. These results had been used to obtain the exact GF employed in [ll[, but the expression of the exact GF is 
valid for all dimensions and it has been used to study FeSi [H, As with the Feynman diagrams, one can rearrange 
for [/ — ?► oo all the diagrams that contribute to the exact Ql^ {k., Zn), by defining an effective cumulant M2^J {zn) that 
is given by all the diagrams of ^^•''(k, Zn) that can not be separated by cutting a single edge (usually called "proper" 
or "irreducible" diagrams). The exact GF ^^•^(k, Zn) is then given by the family of diagrams in figure [T^, but with 
the effective cumulant M2{/(k,z„) in place of the bare cumulant M2^^(z„) = — iD°/(z„ — £/) at all the filled vertices. 

For finite U one has sixteen exact GF 5;{^/(k, z„) that define a 4 x 4 matrix, but when the Hamiltonian commutes 
with the spin it can be split into two independent 2x2 matrices, one for each spin component. We name G-'''^(k, Zn) 
these two matrices, one for each spin component ct, and when there is no possibility of confusion we do not write 
explicitly the a. In Appendix [P] we obtain the expression of G'^-^(k, Zn) as a function of an effective cumulant matrix 
{M} , = M^^((k,z,u) (cf. [D28l) . Similar results are obtained in real space for Qaa'H^^n.) as a function of an 

effective cumulant matrix {M} ^ — M^-jJ {j, z,u) (cf. ID52I) . As before, these effective cumulants define in each case 
two independent 2x2 matrices M 

The derivation given in Appendix [D] is rather general and can be extended to any number of transitions a, but in 
this work we shall only consider the case of [/ — > oo (only one transition per spin a — (0, cr)) and the case of finite U 
(two transitions per spin a = (0, cr), (a, d)). 

IV. THE EXACT GREEN'S FUNCTIONS FROM THE CUMULANT EXPANSION 

1. The exact formal Green's functions for the PAM 

The contribution to ^/^^/(k, z) (cf. Eq. ([59)) ) of the term in the series that has n + 1 effective cumulants is given 
by (cf. Eq. (|D30]) in Appendix ID]) 

{(M . W)" . M}„^, = {M. (W . M)"}_, , (40) 



But note that the usual meaning of vertices and edges is exchanged with that employed in the cumulant expansion. 
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with 

{W}„,„ = W^.^ (k, a, z) = V* (a', k, a) V {a, k, a) G^, (k, z) , (41) 
where we introduced the conduction electron free GF (cf . footnote [S] inAppendix |F]) 

and we have used V{a' , k, a,+) = V (a', k, a) and V{a, k, a, — ) = V^* (a, k, a) (cf. Egs dTBl) and ([27])). To abbreviate 
we define 

A^^, (k, a, z) = (W • M)„„, = ^ (k, a, z) (k, a, z) , (43) 

and introduce it in the series for the exact GF: 

g^^, = |m + M • W • M + M- (W • M)^ + . . .} 

I- J eta' 

= |m + M- A + M-(A)V...| ^ = {M-(I + A + A^ + ...)}^^, . (44) 



We now use 



so that 



(I + A + A2 + ...) =I + A-S=(I-A)"\ (45) 



gii,{k,a,z) = {M-{I-A)-'\ (46) 

which are the components of the matrix G-'^-^ (k, a, z): 

G-'^^ = M-(I-A)"^ =M-(I- W-M)"\ (47) 
We can also express M (k, a, z) as a function of G-^-'' (k, ct, z): we use Eq. (|47| to write 

M = G-''-''-(I-A) =G^-''-G-''-'^-W-M, (48) 

so that 

(l + G^^ • w) -M = G^^, (49) 

and therefore 

M = (l + G-''^ • w) ^G-'^^. (50) 



The calculation of the effective cumulant M is as difficult as that of the exact GF G-''-'^, and to obtain an approximate 
GF we shall substitute the exact M by one that corresponds to an exactly soluble model that is similar to the system 
of interest. To this purpose we shall employ the same Anderson model, but for a conduction band that has zero width 
as well as local hybridization (i.e. e (k, a) — Eq (c) and V (a', 'k,a) = V (a', a)) . This model has the same cumulant 
graphs that appear in the system of interest, but its GF G-f^'°'*{z) can be calculated exactly. To find the approximate 
effective cumulant M"^' (z), we then substitute G-''-^ (k, cr, z) in Eq. (15(11) by G-^ (z) and obtain 

M^P (z) = (l + G^f'""' (z) ■ W^*) • G^f^"' (z) . (51) 

This is independent of k because both G-^ ^'"-^ (z) and 

V*{a',a) V{a,a) 



{W'^*(-)}„,„, 



Z- So (cr) 



are also independent of k. The approximate GF is then obtained by substituting M in Eq. (|47p by (z) from Eq. 

m 



a) 



b). 



+ 



K C7[^ai ki 



aj ^ a'i' 



a 



ttiJi^ a'j' 



— -{y — <y 



a2j2 "'Ji 



k'o-'^o-aka 



^ — xj- — ^ 



"iji^ a'j' 



Q!2j2' 



k'lf^i^^ k 
<i— — f J o- 



= <3- 



gj 



Ail 



k'lfJl^f^ k 
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FIG. 1: Cumulant expansion diagrams of the PAM, that give the one-electron GF in the Chain Approximation (CHA) 
The filled circles (vertices) corresponds to the /-electron cumulants and the empty ones to those of the c-electrons. The lines 
(edges) joining two vertices represent the perturbation (hybridization) (cf. Rules 3.1 and 3.2 in reference 3|). a) Diagrams for 
the /-electron GF G^^^.^, ^,{iui„ = Zn) in the CHA, represented by the filled square to the left, b) Diagrams for the c electron 
GF cri^") in the CHA, represented by the square symbol to the left, c) Diagrams for the f-c electron GF Gj{k,<T( 

Zn ) m 

the CHA, represented by the square symbol to the left. 



2. The exact formal Green's functions for the impurity Anderson model 

In the impurity case, only the Fourier transform of the imaginary time is necessary (cf. Eq. psp ). and to abbreviate 
{{y {f;j,a,u^ -,uj) r (/;j',a',u',w'))+) = g^^ii,a,u = -, cj; j', a', u', w') we define ^^^,(j,iw). Instead of Eq. 
(l39t we then have: 



gff{j,a,u = -,Lo;j\a',u',iu')^gii,{j,t,o) A{u + u') A{u + u') <5(j,j'). 



(52) 



and most of the derivation employed in the previous section can be extended to this case. 

The contribution to Q^^^, (jj, z) of the term in the series that has n+1 effective cumulants is given again by (cf. Eq. 
(|D54p in Appendix |D]) 

{(M . W)" . M}„„, = {M. (W . M)"}„„, , 

where M is defined in Eq. (|D52p : 

{M}^^, =M^{f(j„z„,«), (53) 

and W in Eqs. (|D53ID50p 



As before we define a matrix 



A = W M 



(54) 



(55) 



with components 



A^^, (k,CT,z„) = (W-M) 



(56) 
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and then we obtain a matrix with components Gaa'ih^ ^n)'- 

G^f ^M-il- Ay^ . (57) 

As before 

M = (l + G-''^ • w) ^G-'^-'^. (58) 



and if we substitute G-^-'' (ji, z„) by G^^'°'*{ii, Zn) = G^^'""* (z„) we obtain an approximate effective cumulant for the 
impurity problem. 

M«P (z„) = f I + G^^''^* (2„) • W)"' • Gff''^' (z„) (59) 



Introducing M'^^ (z„) into Eqs. (|55l57p we obtain the approximate GF we shah use in our calculations. 

V. CALCULATION OF THE ATOMIC GREEN'S FUNCTIONS 

When the conduction band has zero width and the hybridization is local (i.e. independent of k), the eigenvalue 
problem of the Hamiltonian introduced in Eqs. (|7l8p has an exact solution [13], and the GFs can be calculated 
analytically. Taking E'k.CT = and VjQ,k,CT — ^q,ct the problem becomes fully local, and one can use the Wannier 
representation for the creation and annihilation operators Cj ^ and Cj.a of the c-electrons. We then write Hr — X]j 
where Hj is the local Hamiltonian 



cr a a, IT ) 



(60) 



and the subindex j can be dropped because we assume a uniform system. 

We shall denote with |n, r) the eigenstates of the Hamiltonian Hj with eigenvalues En^r, where n is the total number 
of electrons in that state, and r characterizes the different states. These eigenstates satisfy 

T-L |n, r) = en,r \n,r) , (61) 

where T-L corresponds to that in Eq. (fTTj) but for a single site, and £n,r = En^r ~ (cf. Eq. (fT2)) ). The states |n,r) 
are usually obtained by diagonalization, and employing Eqs. (|15l6ip we find 

exp[-/3-H]f(7,T)y(7',0)|n,r) ^ cxp (r - /?)] \n" ,r") {n" ,r"\ exp [--Hr] r') (n', r'| Fy |n,r) 

n'n"r'r" 



n n r r 



(62) 



Employing Eq. (jC13p of Appendix ([Cj we calculate the Fourier transform of {^yXj,a{T)x'^- ^,{T'yj 

(X,, „(..,) X]^))^^ = /^{u^s+uj's) ^(x,-„(r)4„,(0))^^^exp[zT u^] (63) 



and we shall then abbreviate 



l^(x,A^s) Xl^,{j,))^^ = A(^. +0^1) GLt^^'i^^s). (64) 
Employing the grand canonical potential Vt — — fcTln^ exp(— /3e„_r) |13il we find for < r < /3 



(x,,„(T)4„,(0))^^^=cxp (/?!]) |g(n,r|exph/3H] X,-„(r) 4„,(0) |n,r)| 



(65) 
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and from Eq. (j62 



(x,-„(r)xj_„,(0))^\ -exp(/3f]) <^ ^ exp [-/3e„,, + (£„,,- e„,,,0 r] (n,ar| X,- „ |n', r') (n', r'| |n,r) 

(66) 

Integrating Eq. (|66l) and using exp [z Ws/3] = —1 and the properties of the Xj ^ we obtain from Eqs. (|63II64I) 



E 



exp(-/3£„_i^r) + exp(-/3g,, 

lUJg -\- Sfi—\ J- ^n.r' 



(n-l,r| Xj^^\n,r'){n,r'\ |n-l,r). (67) 



Equivalent expressions for tj^^'f*(iws), Oj^'^i'^'^iiuJs), Q'^J^f^ {i^jJ s) are easily obtained. One can also consider these func- 
tions as matrix elements of four matrices G^^'°'^{iujs), G'^'^'°'*{iuJs), G-f'^''^*'{iujs) and G'^'^'''*(iws), and one can also 
define a larger matrix G°'^{iLUs) that includes these four matrices, but it is not yet clear whether a formulation that 
uses this larger matrix would have any advantage. We then define 



G'^'ilLJ,) = 



Gff-^^iws) Gf^^^^iujs) 



(68) 



which would be a 6 x 6 matrix for finite U and a 4 x 4 matrix for infinite U . 



VI. DETAILED CALCULATION OF THE APPROXIMATE GF 

A. Introduction 

The hybridization constant V,-,k,CT in Eq. ([2]) is given by Eq. ([3]), 

Vj^i^,a = ^i=K(k)exp(ik.Rj), 

and when the Hubbard operators are introduced into Eq. ([2]) , the hybridization Hamiltonian Hh is transformed into 
Eq. dlD: 



jba.ka 

with hybridization constant VjaMu- The a = (6, a) describes the transition | a >— >■! b >, where the local state | a > 
has one electron more than the state \h >. To simplify the calculation one defines the parameters w(j, a, k, cr, u) in 
Eq. (HH) 

^;(j, a, k,cr, +) = l/,\ba,k,o-, 
i;(j,a,k,cr,-) = V^jtf>a,k,a: 



where m = ±, and in the PAM we use V{a, k, cr, w), defined in Eq. (|27p 

w(j, a, k, (T, 7i) = V{a, k, cr, w)A^s ^ exp(mk.Rj). 

After applying the rules for calculating the GF, it is convenient to return to the explicit use of complex conjugates, 
and we introduce F (a, k, cr) in Eq. (jD16|) 

y(a,k,CT) = l^(Q!,k,cr,u = -), (69) 

so that Eq. (|DT7l) 

F*(a,k,cr) = F(a,k,CT,w = +) (70) 

follows from Eq. (gT]). 

There are four local states | >, | + >, | — > and | d >=| +, — > per site, and there are only four X operators 
that destroy one local electron at a given site. We use the index — 1, 2, 3, 4 to characterize these X operators: 



15 



4 


1 


2 


3 


4 


a — (6, a) 


(0,+) 


(0,-) 


i-,d) 





(71) 



so that /j. = 1,3 destroy one electron with spin up and /j. = 2,4 destroy one electron with spin down. We use 
a — + and a = — instead of a =t and a —\. to emphasize that the spin belongs to a local electron. 
The matrix W, employed in the PAM calculation, is defined in Eq. (jD29p 

{W}^, ^ =M^a',a (k,a,z„), 

and its matrix elements are defined in Eq. (|D26p 

Wa'.a (k, cr, 2:„) = V*{a',k,a)V{a,'k,a) Cj"^ (k, z„) , 

where 2„ = iujn are the Matsubara frequencies. A related matrix appears in the impurity case in Eq. (jD53|) 

with matrix elements defined in Eq. (|D50p 

* k 

The hybridization is spin independent in the Anderson model, so we have 

V{Oa, k, a) = V{ad, k, a) = F(Oa-, k, cr) = V{ad, k, cr) = 0. 

We shall assume a purely local mixing, so that Kr(k) in Eq. ([8]) is k independent, and when we introduce the 
Hubbard operators we obtain 

y(Oa,k,cr) = (72) 

V{ad,k,a) ^ aV, (73) 

where we have also assumed that Vcr(k) is independent of cr = ±1. 

As discussed in the introduction of Section IllH when the Hamiltonian is spin independent or commutes with the z 
component of the spin, the 4x4 matrices G-^-'', M, W and A = W.M can be diagonalized into two 2x2 matrices, 
e.g.: 

(74) 

1, 3 appear in G'^'^ and Ix =2,4 

(75) 
(76) 

(77) 
(78) 



In this matrix the indices defined in Eq. ([7T|) have been rearranged, so that 
appear in G| . 

Employing Eqs. ([7^r75|) we find for the PAM (cf. Eq. (gl])) 

w^(k,z)= \v\' gl^{k,z) 
W,(k,z)= \v\' Gl^iKz) (^^^ 

For an impurity located at the origin we find 

w^{z)^\vu^{z) r M, 
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where 



W = ^c%(k,z). (79) 



k 



For a rectangular band with half- width D in the interval -B], with B = A + 2D we find (the minus sign of 
Gc,a- (k, z) is included in the logarithm) 

, ^ 1 , z — B + u , , 

2D z + A + fi 

where the /i appears in (pa{z) because of the e (k, a) = i?k,cr — M in ^/c a C*; •2) ■ 

Both for the PAM (Eq. (|47)) ) and for the SIAM (Eq. (ISTl) ) we have the same relation for the submatrices in Eq. 
dZll): 

G^-^ = M,. (I - A,)-\ (81) 

and as before (cf. Section HIT)) 

M,= (i + G{-^.W,)"'.G//. (82) 



B. The approximate G^^'"*' GF for the Periodic Anderson Model 

The calculation of the exact effective cumulants Mg. is as difficult as that of the exact G^-'^, and the atomic approach 
consists in using instead the M^- of a similar model that is exactly soluble. Taking a conduction band of zero width 
aX £0 = Eq — ^ and a local hybridization, the PAM becomes a collection of independent atomic Anderson systems 
that can be solved exactly, and we call this solution the atomic Anderson solution (A AS). We then employ the A AS 
to calculate the corresponding exact Green's function G^'^'''*(z), and we define an approximate effective cumulant 
by introducing this GF into Eq. ((82|) using (zo;) = —1/ (iw — £0) in the corresponding W°*. This procedure 
overestimates the contribution of the c electrons because we concentrate them at a single energy level Eq^ 

and to moderate this effect we replace V'^ by in the calculation of the M"'', where A = tt /2D is the Anderson 
parameter (but we keep the full \ V\ when we substitute this into Eq. (IHlT) to calculate the approximate G^-^^^^^). 

1. Imaginary frequency and reciprocal space 

To calculate the approximate G^-'^'°^(k, iw) for imaginary frequency and reciprocal space we then introduce the 
into Eq. dH]), with A^p= W^.M^^' and W,, defined in Eqs. ([75T761) with the fuh We can now write 



mil 

,TO3i m33 , 



m22 

, W42 m44 , 



(83) 



and 



On ai3 

, 031 033 



022 024 
, 042 033 



(84) 



and from Eqs. (|75I76I) we find 



On tJ"^ (k, z) (mil + m^i) 
033 = (k, z) (TO33 + 77113) 
031 = ail 
ai3 = 033 



022 = 5°; (k, z) (m22 - 77142) 
044 = 5°; (k, z) (m44 - TO24) 
042 = —022 



(85) 



024 



-044 
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To derive our approximate G^-^-^p we substitute both and A°p= W^-.M^p into Eq. dH]), and we obtain 

/ mil TOi3 



and 



(k, z) (miimgg - migmgi) 



G(^'''^(k,itj) = ^ — ^ (86) 

1 - \V\ g^.^ (k, z) (mil + + TO13 + TO31) 



I 77122 "^24 \ |^|2^0 (-J^^^) (^^22^^^ _ ^^^^^^-j / 1 Ij 

Gf /-^(k, ..) ^ V"^" "'"^ , „ ^ V^. (87) 

1 - g"^i (k, z) (to22 + ^44 - TO24 - TO42) 

Note that in this approach the 'M.^P{iuj) are independent of k, and that all the k dependence of our approximate 
Glf'''P{k,iuj) comes through the Gl^ „ (k, iw) in the Eqs. (|86ll87p 



2. Real space and imaginary frequency 



In the previous section we derived the reciprocal space and imaginary frequency GF for the PAM in the atomic 
approximation, but sometimes it is necessary to use the GF in real space, with the / electron being created and 
destroyed at the same site. These GF are denoted by G{,^ {iuj) and are given by 



k k 

Considering a rectangular conduction band and transforming the sum into an integral, we obtain 

-D-^l 

-D-fi 



(89) 



where ^ is the chemical potential and D is the half-width of the conduction band. Substituting Eqs. ([86l and [87]) into 
Eq. ([89]) we obtain that 



and 



\'-^31 '-^33 / 



M°^- |y|^t/°>.(£,iw)(TOiiTO33-TOi3m3i) ( ^ } 

1 r+JJ-t^ ' \ — 1 1 

— / de 2 ^ 

2D ,l_n-,, 1 - tj">.(e, iw) (mil + 77133 + mi3 + msi) 



(90) 



Gl^iiuj) 



^22 "-^24 



1 1 



1 

2D 



M'^P-lVl ^?°^(e, iw) (m22m44 - m,24m42) ^ 

1 - l^^l^ Scii^, i^) (?7l22 + m,44 - m24 - m42) 



de- 



The variable £ is present only in Q^i ^ and the integration is straightforward. We find 



(91) 



Gf^(z) = M^^' 



z- D + a + V^Mi^ \ 
2D \ z + D + fi + V^M(f J [ ^ ^ 



1 -1 
-1 1 , 



(92) 



where 



G{^(z) =M7 + :^/n 



2D~ 



'z-D + fi + V^Ml^\ 



D 



= mil + mi3 + msi + m33. 



1 1 



e. 



(93) 



(94) 



'I 

and 
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M/^ = TO22 + - m24 - (95) 



9-t- = miirnaa - misWai, (96) 



0; = "l22m33 - m24TO42- (97) 

In Appendix [E] we define and calculate the formal expressions of the matrices G^=^°P(k, iw) (cf. Eq. (|E5IE6|) ) 
and G^^''^P{k,iuj) (cf. Eq. (jElllE12|) ') associated to the crossed GF of the impurity, as well as the GF of the pure 
conduction electron GJf '"^(k, iw) (cf. Eq. (jE16IE17[) ). We can also describe the conduction electrons in real space 
and imaginary time: the corresponding G^'^''^^ (iuj) are given in Eqs. (|E18|E20|) and the G'^^''^P{iu}) in Eqs. (|E21IE22p . 
In a similar way we obtain G^'^'°-P{iuj) (cf. Eq. (jE23| ). and we can use this relation to express all the other GF (cf. 
Eqs. (IE24IIE27IE28IE29I) ). 



C. The approximate G^^ "'' GF for the Impurity Anderson Model 

As in the case of the PAM, we substitute the G^^ in Eq. (|82p by the exact solution G^-'^'"* of the problem with 
zero band width and local hybridization, and obtain the corresponding effective cumulant M"^'. The conduction band 
corresponding to this approximation then has i?k,cr = Eq, so that for all values of k it is (k, z) ~ — 1/ (z — Eq) , 
where So = — ^. From Eq. ([75]) we then find faiz) fai^) — ~ ^o), and substituting into Eqs. (|77l78p 

we obtain the W^p that should be used in Eq. ^ to calculate M^^. To define the approximate GF Gl^'^^Piiuj) 
introduced in the present work, we substitute this approximate MJ^p into the exact expression Eq. (I8ip . but in this 
equation we use the that corresponds to the conduction band with full width. The Eqs. (|83|84|) for M"^ and 
A^P= Wo-.M"'' are also valid for the SIAM, but using Eqs. (|77l78p we find different expression for the Uij: 

an = ip-f{z) {mil + msi) ; 022 = </?i(2) ("122 - "^42) 

0-33 = Vti^) + mis) ; 044 = </?i(z) (^44 - "^24) 

^31 = On ; 042 — —0,12 

oi3 = 0,33 ; 024 — —044 

If we now substitute these M"^' and A°p into Eq. ([8T|) . we obtain the approximate G^^'^^p for the SIAM: 

/ mil 7^13 1 |T^|2 \\ f \ 

\-\V\ tp-f((iu}))[miim33-mi3m3i) 

ffn„, . \m3i m33 J . , , ^ 

G^'^Piiuj) = ^ ^ 2 ^ (99) 

1 - |y| ip^{{iuj)) {mil + m33 + 17113 + m3i) 

and 

/ TO22 m24:\ |Tr|2 //. NW 

- 1^1 ¥';((*W)) (TO22m44 - TO24m42 

Gf -(^0.) = ""^ , (100) 

1 ^ 1^1 <y5i((«^^)) (^7122 + m44 - TO24 - m42) 

In Appendix |E] we define and calculate the formal expressions of the matrices Gl'''''P{j, ^ 0,k,icj) (cf. Eqs. 
(|E42IE43p ) and G''/'°'P{k,y = 0,iuj) (cf. Eqs. (|E48IE49p ) associat ed to the crossed GF of the impurity, as well as 
the GF of the pure conduction electron GJ^'^'"''(k, k', iw) (cf. Eqs. (|E53IE54p ). We can also describe the conduction 
electrons in the Wannier representation: the corresponding Gl'^'°-P{iuj) arc given in Eqs. r??lE57l) and the G^J-''P{iuj) 
in Eqs. (jE58IE59|) . In a similar way we obtain G^'^'°-p (iuj) (cf. Eq. (jE60[) ). and we can use this relation to express all 
the other GF (cf. Eqs. (IE6mE64IE67lE68l) '). 





Appendix A: The sign of the contribution of a graph. 



For convenience we summarize here the Appendix C of reference Q, where one should hav^ v{a, k, a, ±m) instead 
of V{a, fc, cr, ±u) in the item 2(b) of Rule C.2. 
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Here we discuss the sign that must be given to the contribution of a given graph, and we are only interested in the 
case without external fields, i.e. with ^ = 0, but we shall keep some results for ^ 7^ that are convenient to understand 
^ = 0. The rules for drawing the graphs that appear in the calculation of the averages • • • Y{ln))^)^ are given 

in Rules 3.3 and 3.4 in Q. In item 4 of those rules, the Fermi type lines running to each vertex were paired in an 
arbitrary way for ^ = 0, and several open and closed loops were formed in this way, where all the open loops must 
have two external vertices. A definite sense was arbitrarily assigned to each of the loops, and we call this direction 
the "sense of the loop". In the following discussion we consider only Fermi-type operators, because the position of 
the Bose-type operators does not affect the sign of the contribution, and we shall mean Fermi-type operator when we 
say "operator" in the remaining of this Appendix. It is now convenient to introduce two concepts that shall be useful 
in the present computation. 

Definition C.l. 

A graph is in a "perfect ordering" when the following relations are satisfied: 

1. For all the open loops, r increases in all the vertices of the loop in the sense of the loop. 

2. For every closed loop, r increases in the sense of the loop for all the vertices but one (it is impossible to satisfy 
(1) for a closed loop). 

3. All the r in a given loop are either smaller or greater than all the r in all the other loops of the graph.* 

There are many ways to choose a perfect ordering of a graph, but the particular choice is not important provided 
that we use always the same one after it has been chosen. 
Definition C.2. 

Several Fermi-type operators of a graph contribution are in a "perfect order" when: 

1. The y-operators are written from right to left following the perfect ordering we have chosen for their graph. 

2. For the two operators of each internal edge (they have the same r) we write the X-operator to the left of the 
C-operator.» 

As a starting point we shall consider rules that are also valid for ^ 7^ because they are simpler to state although 
less systematic to apply than those given by Hubbard Q for ^ = 0. We shall consider explicitly the graphs for 

(^(y{1) ■ ■ ■ Y{r)^ ^ , i.e. a GF with r external operators Y{1) ■ ■ ■ Y{r) (but the rules are also valid for vacuum 

graphs). The n-th order term of this GF contains the average 

^ (y (1) • • • Y{r) (H'in) ■ ■ ■ i?'(r„))) (Al) 

and the application of Theorem 3.1 from f§\ to this equation gives all the 71-th order graphs^. 
Rule C.l. 

To obtain the sign associated to a given graph, multiply the parities of the following two permutations: 

1. It takes the operators from the order used to write Eq. (|Aip into the perfect order. 

2. It takes the operators from the perfect order into the order in which they appear in the final expression that 
gives the graph contribution. 

As the operators Hh are of the Bose type and can be moved freely inside the ordered parenthesis in Eq. (|Aip . in the 
first step it is necessary to consider only the permutation that takes the external Fermi-type operators to their perfect 
order. The Rule C.l. is just the application of Theorem 3.1 in Q in two steps, and the only reason to proceed in 
this way is that the perfect order of the ^-operators in a graph provides a reference frame to organize the calculation. 

For ^ = we shall give rules with the same labels employed by Hubbard Q , because of their similarity. In this case, 
there is only an even number of lines running into each vertex, and for any CV this number is two .This simplifies 
the treatment, and the first step is the same step 1) employed in Rule C.l: this is just rule "d" of Hubbard. 

To calculate the change of sign that corresponds to step 2) in Rule C.l we proceed in three steps. 

First we consider all the open loops that pass through each vertex, and note that in the perfect order, the X- 
operator is to the left of the C-operator in all the internal edges. To be able to pair operators of the same type at each 
vertex (otherwise the corresponding cumulant vanish) it is necessary to change the order of these two operators (with 



In the absence of an external field (i.e. when ^ = 0) the H'{t) in this equation is equal to Hi^ij) (cf. Eq. 



20 



a change of sign) when the arrow in the edge points towards the CV. To correct for the sign missing in Eq. ([T6)) one 
must also add a factor ±u to the v{j, a, k, a, ±u) in Rule 3.5.2.C in Q, and these two factors correspond to Hubbard's 
rule "b" . In the present problem there are only two edges at each CV, and when both are internal, the effect cancels 
out and the rule is not necessary. To prove this result, note that according to Rule 3.6b' in Q, the cumulant 




at each CV, is already written with the K-operators in the perfect order, with the —U2 corresponding to the outgoing 
arrow. The contribution to Rule 3.5.2.C in [3| of the two internal edges running into the CV after correcting for the 
missing sign in Eq. (|16|) . is then 



{+U2)v{j2,a2,ks,(T2,+U2){-Ui)v{ji,ai,ks,<7i, -Ml) (A3) 

As there is particle conservation, we have ui — U2, and when we multiply into the minus sign due to the exchange 
of the X with the C operators on the line with the arrow towards the CV, the overall sign is always plus. Hubbard' s 
rule "b" is therefore not necessary for all the CV with two internal lines. 

For any CV with only one internal line (and therefore one external line also), one must multiply the w(j, a, k, a, ±u) 
into ±M and also into —1 when the internal edge points toward the CV . This is the only effect that remains in the 
PAM of Hubbard's rule "b" . 

The discussion above fails for closed loops because r can increase in the sense of the loop in all vertices but one. 
After putting all the operators in perfect order and then exchanging the X-operator with the C-operator for all the 
lines with arrows pointing to a CV, the first and last operators in the resulting expression belong to the same vertex, 
and should therefore be brought together . These two operators are separated by an even number of Fermi operators, 
but bringing them together by an even permutation would still leave them against the order of the loop, i.e.: the 
operator at the left would correspond to the edge with the arrow pointing toward the vertex. A permutation of odd 
parity is then necessary to put all the operators of any closed loop in perfect order, and this is Hubbard's rule "c" . 

After the three steps discussed above, the F-operators that were in the order given by Eq. (jAll) are now paired 
at each vertex according to the loops of the graph considered, each pair written in the sense of the loop. We shall 
denote with (as,/3s) the two indices of the F-operators of each of those pairs, written already in the sense of the 
loop, i.e.: /3s — Q^s- All the pairs that correspond to a given vertex are still separated by many pairs that belong to 
other vertices of the graph, but it is only necessary an even permutation to put together all the pairs of each vertex. 
The pair associated to each CV is already in the same order of the cumulant of Rule 3.6 6' in Q, and only remains to 
consider the cumulants associated to the FV. If there are p loops crossing an FV, we have already the corresponding 
operators in the order (ai,/3i), • • • (ap,/?p) while in the cumulant associated to that vertex by Rule 3.6. (2). a in [3| 
they are written in the order 1^(71) • • • ^(72^), where 71 • • • 72p, correspond to the same(ai, /31), • • • (ap, /3p) but in a 
different order. It is then necessary to associate to each of these cumulants a ± given by the parity of the permutation 
that takes (ai, /3i), • • • {up, (ip) into 71 • • • 72^. This is Hubbard's rule "a" . 

It is now convenient to put together the rules for the calculation of the sign required by Rule 3.7.(2).d or Rule 
3.7a.(2).e. 

Rule C.2 

To calculate the sign of a graph with f = 

1. Define a perfect ordering for the graph according to Definition C.l. 

2. The sign of the graph is the product of the following factors 

(a) When there are p loops crossing an FV, denote with (as,/3s) the indices of the two X-operators of the s-th 
loop at that vertex (s — 1, • • • ,p), written already in the sense of the loop (i.e.: (3s a^)- The 2p Fermi- type 
operators at that FV appear in the cumulant of Rule 3.6a in the order 1^(71) • • • y(72p), where the 71 • • • 72p are 
the same (ai,/3i), • • • (ap,/3p) in a different order. For each FV multiply into a ±1 given by the parity of the 
permutation that takes(Q!i, /3i), • • • (a^, /3p) into 71, • • • 72p- 

(b) For any CV with only one internal edge multiply the V{a, k, a, ±u) of Rule 3.6. (2). c into (±u), and also into a 
further — 1 when the arrow of the internal edge points toward the C V. 

(c) There is a factor —1 for every closed loop. 

(d) If the graph is employed to calculate a GF with r Fermi-type operators written in the order Y{1) ■ ■ - Yir), 
multiply into a sign given by the parity of the permutation that takes (Y (1) ■ ■ ■ Y (r)) into the same operators 
written in the perfect ordering chosen for the graph. This item does not apply to vacuum graphs. 
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Appendix B: Counting graphs and the symmetry factor. 

For convenience we reproduce here, with very minor changes, the Appendix D of reference 

As discussed in appendix |X] of reference [3|, the n-th order term of the perturbative expansion of the GF 
Y{1) ■ ■ ■ Y{r) ] j contains the expression in Eq. (lAip of [3], and its contributions have the form 



H 



(-1)' 



V{lnAn){{Y{h)Y{l',)---Y{l^)Y{Q)^f (Bl) 

(cf. Eq. (3.11) in but with the r external operators Y{1) ■ ■ ■ Y{r) inchidcd in the averages. When the Theorem 
3.1 in [3] is apphed to these averages, the n-th order contribution can be associated to a family of graphs, and many 
of them are disconnected and composed of several connected graphs. We label each topologically distinct connected 
graph with an index a, and we use tIq, to denote the number of times that the a graph appears in the nth-order 
graph. It is clear that there might be several identical contributions associated to the same n-th order graph because 
all the n\ permutations of the edges of a given graph give the same contribution. These identical contributions should 
be counted as different contributions every time they correspond to a different partition in cumulants. The correct 
number of times that a topologically distinct graph of n-th order gives the same contribution is then 

n (B2) 

where is the symmetry factor of the connected graph a and is calculated using the Rule D.l discussed below. To 
derive this result one applies the same arguments employed in Ref. the factor 2" of that reference is not present 
in our expression because the pair of vertices of any internal edge can not be exchanged (cf. the definition of the 
coefficients of Eq. discussed after Eq. ([TO| ). 

To calculate the symmetry factor ga it is enough to adapt the rule given by Hubbard in Ref. 0, Appendix B 
. The calculation seems rather obvious in simple cases, but it is convenient to give the rule to deal with the more 
complicated ones. 

Definition D.l 

A vertex is said to be "internal" when all the lines running to it are internal lines.* 

In the PAM, only Fermi lines can run into an internal vertex, because of the form of the interaction (cf. Eq. ©). 
Rule D.l 

To calculate the symmetry factor g of a connected graph with p f and vertices FV and CV respectively: 

1. Number the FV with 1, 2, • • • ,pf and the CV with 1,2,... ,pc so that 1, 2, • • • ,qf correspond to all the internal 
FV and 1, 2, • • • , to all the internal CV. 

2. Form the pf x pc matrix N, with elements Ni j, where Ni j is the number of Fermi edges joining the FV i to the 
GVj . 

3. Let gi be the order of the group of permutations Vi of the qf x Qc ordered pairs («, j), which has the property 
that if any permutation of Vi is applied to the indices i — 1,2, ■ ■ ■ ,qf and j = 1,2, - ■ ■ ,qc oi the matrix N , this 
matrix is left unchanged. 

4. The symmetry factor is then 

9 - 9lf[f[iN^.,l) (B3) 

Appendix C: The Fourier transform of Green's functions in imaginary time 

We are interested in the GFs defined in Eq. with only two operators Y{'y, r): 

e(7i,n;72,T2) = ^(f(7i,ri)y(72,T2))^^^, (CI) 



7,1 
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as well as in their Fourier transforms. Introducing 

f} = -iln[exp(-/3H)], (C2) 

we can write for ti > T2 (cf. Eq. (fTSl) ) 

^^(71, Ti;72, T2) = exp (jSD,) Tr {exp {-13%) exp {t{H)Y^^ exp (-r^'H) exp {t^T-L)Y^^ exp (-r2'H)} (C3) 
and using the properties of the trace we obtain: 

^(71, n; 72, T2) - exp ipn) Tr {exp exp ( (n - T2) 'H)Y^, exp (- (n - T2) } = 

= ^(f(7i,ri-r2)r(72,0))^^ ^ F {t^ - T2) . (C4) 

Similarly for ri < T2, we have (because both are of the Fermi type) 

^(71, n; 72, T2) = exp (/3ri)Tr {exp (-1)^(72, r2)y(7i, n)} = 

= — exp (^f2) Tr {exp {—(3'H) exp (t2'H)F^2 exp {—t^'H) exp (r^'H)Kyi exp (— t^H)} = 

= - exp Tr {exp ( (n - T2) H)r^, exp (- (n - ra) exp F^J = 

= - exp ipn) Tr {exp (-/3H) exp ( (n - T2 + /?) H)^^, exp (- (n - T2 + /?) H)^^, } , (C5) 

and as Ti — T2 + /3 > we have 

e(7i,Ti;72,T2) = -T(Ti-T2 + /3), (C6) 

and finally 

g(7i, n; 72, T2) = (ri - r2) T (n - T2) - (r2 - n) F (n - T2 + /3) . (C7) 
The time Fourier coefficients are given by 

C/(7i,wi;72,i:J2) =^y" dri j dT2 exp (cjiTi + ^2X2)] ^ (^F(7i, Ti)r(72, T2)^ = 
1 /-^ /-^ 

= - / dTi / dT2 exp [i (wiTi + ^2X2)] (6* (ti - r2) F (ti - T2) - 9 {t^ - (n - T2 + /?)) . 
P Jo Jo 

(C8) 



We change variables 



and we find 



X = Ti - T2 

y = TI (C9) 



1 



/3 /-/J-y 



^(71, wi;72,a;2) = - y dy exp [iy + 1:^2)] / dxQyj^{ixLo^\{Q (x) F {x) ~ Q {-x) F {x ^ fi)) = 



y 

— j dyexp [iy (cji + 0^2)] ■( / dxF (a;) exp [ixwi] — / dxF [x + f5)eyi-p[ixu)i\\ . (CIO) 



1 



/3 7o 

Changing variables x + /3 = 1^ in the second integral we obtain 



/ dxF {x + (3) exp [ixwi] — I 

J-v J 13 



(0 exp [ixcji - ;9)] (Cll) 



and employing exp [ia:a;i/3] = — 1 we hav^ 

f/3 r /-/S-y 



^/(7i,'^i;72,W2) = ^ y dyexp [iy (wi + a;2)] |y da;i^ (x) exp [ixwi] + y dxF (a;) exp [ia;a;i] | 

= A(wi+ti;2) / dxi^ (x) exp [zxwi] (C12) 
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where loi and 0J2 are given by Eq. (|19l) for Fermi-like operators. When the two operators 1^(7, t) are X-operators we 
can write 



= A(L.i+^2)^''rfri^(x,,,„,(ri)xj^„^(0))^^ exp[irit^i] (C13) 

Appendix D: The exact GF as a function of effective cumulants. 

1. Rules for reciprocal space and imaginary frequencies (Valid for the PAM) 

In Section [IV II we use the exact solution of the model with zero band- width and local hybridization to approximate 
the effective cumulants that appear in the exact GF. Here we shall use the prescriptions given in Rule 3.7 of Section 
|TT]to calculate the diagrams in imaginary frequency and reciprocal space, to obtain the formal expression of the exact 
GF 

^ k, a, u = uj,)Y{f- k', a', u', uj'^)) ^ = 0^, (k, tu;,) A{u + u') A (^.k + u'\^) A {uo, + ^) (Dl) 

in terms of the corresponding effective cumulants, where ujj and ui'j are the Matsubara frequencies (cf. Eq. ll9p . As with 
the Feynman diagrams, one can rearrange all the diagrams that contribute to the exact Q^J^,{'k,iujj) by introducing 
effective cumulants M'^^J, (k, ujj ) , defined by the contributions of all the diagrams of QI/^, (k, iujj ) that can not be 
separated by cutting a single edge (usually called "proper" or "irreducible" diagrams). The exact GF Qlf^,{k.,iujj) is 
then given by the family of diagrams in reciprocal space corresponding to those given in figure [1^ for real space, but 
with the effective cumulants M^^f (k, Ljj) replacing the bare cumulant Af^^,(k, Wj) = —6aa' Da / (iujj ~\- £a) at all the 
filled vertices (here Sa — £b — £a when a — lb, a), cf. appendix IF)) . To calculate the contribution of the diagram with 
n + 1 effective cumulants we follow the steps in Rule 3.7: 

1. We label all the diagrams containing n + 1 effective cumulants M^l^^ (ki, wi, ui; k2, uj2, "2). Because of Eq. ([22]) . 
this effective cumulant is proportional to A (wi -I- W2), and because of the 5 in Eq. (|28|) it is also proportional 
to A (uiki -|- U2k2). The particle conservation also requires a A (ui -f- U2), and the labels we use are shown in 
figure [2] 

We make the product of the following factors 

(a) All the A (u + u') A [uk + u'k') A (oj + lo') that appear in Eq. (jPlI) remain with the effective cumulants, because 

they correspond to all the proper diagrams of Q^J^,{k.,iuj). We then have 

N, A (u'k' - liiki) A {u' - ui) A {J - uji) M^(^,(ki, -tJi, -m; k', J, u') 
xN, A «k'i - U2k2) A {u[ - U2) A {uj[ - ^2) M^^, (k2, -^2, -U2]K,ij[,u[) 
xiV, A (4k'2 - uakg) A (4 - U3) A (4 - W3) M^{^, (kg, k'2, w^, ^2) 

xiV^ A «_ik'„_i - Wnk™) A «_i - ii„) A -a;„) (k„, -w„, -m„; k^_i, 

xiV, A(uX + ^ik)A« + 7/)A(^;+c^) M:{{(k,a;,w;k;,w;,<). (D2) 

(b) The contribution of the n cumulants of conduction electrons 

1 



-f ui e (ki, (Ji) 



(5(ki , k'l )(5(ui , -u'l )(5(cri , )(5(wi , ) 



5(k„,k^)(5(M„,05(cr„,cr^J(5(w„,w^J. (D3) 



iw„ -I- u„ e (k 
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FIG. 2: Expansion of the exact GF employing effective cumulants. The figure represents the collection of all the diagrams with 
n + I effective cumulants: we write k , a , it , a; = kg, Qfo, ^O) '^0 s-^d k, a, ii, cj — k,i+i, ctn+i, Un+i,Cjn+i 



(c) The contribution of the 2n interaction edges 

V {a[ ,k[,(T[,u[)V {ai,ki,(Ti,-ui) 

). (D4) 

There is a factor (—1) for each interaction parameter F (as, k^, ds, u^), and it cancels out for the (k, iw), 
and the 5^^, (k, iw^), but for t?^-^, (k, iw^), and g^'^,{'k,iujs), one of these factors (—1) remains and a change of 
sign is necessarjl2l. This sign is not necessary in Eq. (|D2p because it cancels out like in the t?;^;^, (k, icj). 

(d) A factor ±1 obtained employing the rules given in Appendix El 

The graphs represented by figure [2] can be considered to be in the perfect order of Rule C.l, and we can apply the 
first step of Rule C.l. without any changes in sign. We assume that the contributions to the effective cumulant 
M^^'f (k, cjs) have been already calculated with their correct sign, and we only have to discuss the diagram with 
n + 1 cumulants in figure [5] as a whole. As discussed in Appendix |^ is is not necessary to introduce any sign 
change for all the CV joined by two internal lines, and only those joined by only one internal line (and therefore 
joined also by an external line) have to be considered. These CV appear only in the three GF of the type 
Q^^,{'k,iuJs), G^aiO<.,iujs), and C/;^^, (k, icj^), and in disagreement with previous results [lB|, no change of sign is 
required for these three type of GF, so that in all cases we have only to multiply into 



+ 1. (D5) 

(e) A factor 1/g calculated from appendix |B] We assume that all the factors that appear in the contribution of the 
diagrams corresponding to the effective cumulant Af^^'f (k, oj^) have been already included in the M^'^J (k,uJs) 
itself, and we then only need to use the g corresponding to a chain, that is g = 1. We can then write 



(D6) 



(f) From the two FV external lines we have an extra factor 



3. We now have to sum the products with respect to: 

(a) the momenta kg, the frequencies ujs and the indices Us of all the internal edges, and also divide each sum over 
momenta into ^/N^. This last contribution cancels exactly n factors Ns,oi all the n + 1 factors Ng that appear 
in Eq. (jD2[) . The extra factor Ng is canceled by the spatial Fourier transform of the external lines, and is taken 
care by the item 2.f of rule 3.7. When there is an external line running to a CV, there is a 1/^/7^ associated 
to the internal line running to that CV, but there is no l/y/Wl associated to the external line, because the 
corresponding operator has been already introduced in k space. 
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Because of the delta functions in Eqs. (ID2ID3P we have 

= Us = u'^ = —u, (D8) 

k'= k, = k; = k, (D9) 

and the Matsubara frequencies 

uj' = uJs = uj's = -uj (DIO) 

for all s = 0, 1, . . . , n. 

(b) We have the sum over all as, a'^,, for s = 1, 2, . . . , n and we shall use matrix notation to simplify the calculation. 

(c) Because of the S{as,(j's) in Eq. (|D3P and the spin conservation in the effective cumulants when [az'jH] — (even 

if a's 7^ tts+i) there are no sums left over the labels as- 

We shall now consider the contribution of the factors in Eq. (|D2P : employing Eqs. (jD8IID10|) we can write for s = n 
M^;^-f (k,iLj = -LOn,u = -M„;kjj,w,'j = iOn,u[^ = Un) = Af^^{(k,iLj,u;k, iLj„ = -io,Un = -u); (Dll) 
for s = 1, 2, . . . , n — 1 

KiUo.'^ {^s+i,-uJs+i, -Us+iX,^'s, <) = KiUo.'^ (k, u- k, -w, -u); (D12) 

and finally 

Ml{l, (ki , -^1, -zii; k', u:', u') = M'J/^, (k, u; k, -cj, -u). 

If we now define 

M^/{(k, tuj, u) = ArJJ (k, UJ, u- k, -UJ, -u) (D13) 

we can write the contribution of the n + 1 factors in Eq. (|D2I) in the following form (cf. 3 (a) for the cancelation 
of n factors Ns ) : 

Mll[ (k, lu, u)Ml[l,^_^ (k, lu:, u)... A'^li, (k, tc, u)M^J/^, (k, it,, u) (D14) 

We still have to include in the contribution of the diagram with n + 1 effective cumulants all the factors from Eqs. 
(|D3|D4p . Employing Eq. ^ and Eq. ^ we have 

V{a,k,a,-u) :^V*{a,k,a,u) (D15) 

and it is then convenient to write 

V{a,k,a) = V{a,k,a,-) (D16) 

V*{a,k,<j) = V{a,k,a,+) (D17) 



As before we use Eqs. (jD8|D9P and the conservation of a to simplify Eq. (|D4|: 

V {a's, k's, a's, u's)V (as, ks, as, -Us) = V {a's,k,a,Us)V {as,k,a,-Us) , (D18) 
and we combine 2 (b),(c) introducing Wa,a' (k, cr, w; iws) 

Wa' a (k, a, u ; iuj) = V (a' , k, a, u') V (a, k, a, —u) ; — r . (D19) 

Employing Eas. (jD16ID17p wc then obtain 

W^.,^ (k, a, u! = +; lUJs) = y* {ol, k, a)F(a, k, a)-. ^— (D20) 

iLkJs + £ (.K, a ) 
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and 

W^,,^ (k,a,u' = ^-iws) = y(a',k,a)F*(a,k,(7). 5— -. (D21) 

iLOs — e (k, crj 

We return to the GF (k, w) defined in Eq. (jPlI) . wliicli corresponds to u' = +, so that the factors from Eqs. 
(|D3ID4p could be put then in the form (cf. the labels in figure [2]) 

(k,CT, +;«CJ„)M^Q^_^,a„_i (k, O-, +;jWn-l).--W^a^,a2 (k, Cr, +; icJ2) W^q'^ ,qi (k, cr, +; i^i) , (D22) 

We now introduce the c-electron free GF C/^^. (k, z„) (cf. Eq. (|42p ) 

and from Eq. (|D10p we have Ws = — w, so that 

1 -1 



Ql„{V,-^UJ:)^QlA^M■ (D24) 



iujs + e (k, cr) —iujs — £ (k, ct) 
We can then substitute in Eq. (ID20I) 

l^a^,,, (k,a,+;iL^,) = T/*(a'„k,a)F(as,k,a) (k, iw) (D25) 

and then, to calculate the contribution in Eq. (ID22p . it is more convenient to use the quantity 

(k, a, z) = k,a)F(a,k, a) 5° , (k, z) (D26) 

where now the complex variable z takes the place of the Matsubara frequency iuj. The contribution in Eq. (jD22p 
takes the form 

Wa^,a„ (k,cr;z) (k, ct; z) . . . (k, cr; z) W"„'^,„^ (k,CT;z). (D27) 

To simplify the calculation we now introduce the two matrices (cf. Section [IV ip 

{M}^^, =M^^J(k,z,u), (D28) 

and 

{W}^,^ = W^„.,„(k,a,z). (D29) 

The contribution of the diagram with n + 1 cumulants is then (using the Einstein convention of sum of repeated 
subindexes) 

(+1) X i X {M}„ {W} {M}„^„, {W}„, . . . {W}„, {M}„^„, {W} , {M}„^„, = 

{(M . W)" . M}„ = {M . (W . M)"}„ (D30) 

2. Rules for real space and imaginary frequencies (Valid for the impurity) 

We shall now obtain the formal expression of the exact GF (^(Y{f;j,a,u — — , Ws)F(/; j', a', u', ^ in terms 

of effective cumulants for the system (SIAM) with a single impurity at site j^, and we shall use the prescriptions 
given in Rule 3.7a of Section to calculate the diagrams in imaginary frequency and real space. The diagrams are 
topologically the same employed for the PAM, and we write 

g£,(j,ZL.„u;j',»c^^,u') = (^(Y{f;j,a,u,iOs)Y{f;y,a',u',iu',))^'^ , (D31) 
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but as there are local / states only at the site j^, we must then have j' = j = jj, and we write 

gii,ij,iu;s,u^ -,y,tuj'^,u') ^gii,{j,,tLOs) Aiu + u') 5(j,J) 6 A (lo. + uj'J. (D32) 

As in the previous Section|DT]one can rearrange all the diagrams that contribute to the exact Gaa' (j, i^^s, j', iuj's, u') 
by introducing effective cumulants M^;^, (j, ioj^, u; j', iw^, m'), defined by the contributions of all the diagrams of 
Q^J^,{],iujs-,u]]' ,iuj'^,u') that can not be separated by cutting a single edge (usually called "proper" or "irreducible" 
diagrams). The exact GF fj^;^, (j, iws, u; j', icj^, u') is then given by the family of diagrams in figure [T^, but with 
effective cumulants M^^^^ (ji, wi, -ui; j2, 1^2, "2) replacing the bare GF Q^j ^^,{}si^^s) at all the filled vertices (as usual 
£a — £b — £a when a = {b,a), cf. appendix |F|. To calculate the contribution of the diagram with n + 1 effective 
cumulants we follow the steps in Rule 3.7a of Section Hi Bl 

1. We label all the diagrams that appear in the expansion of the ^^^^^ (ji, wi, ui; j2, a;2, 1(2) corresponding to the 
SIAM, and in figure[3]we show the diagram that contains just n+1 effective cumulants Af^^£^ (ji, wi, ui; j2, W2, ^2). 
Because of Eq. ([22|) . this effective cumulant is proportional to A {uji + UJ2), the particle conservation requires a 
A (ui + U2) and for the case of an impurity at site the contribution is also proportional to S (ji, ji) 5 (j2, 
because there are / states only at that site. We shall then use the GF G;^^, (ji, iujg) defined in Eq. (jD32[) . 

2. We make the product of the following factors 

(a) All the A{u + u') 5 (ji, j) 5 (ji, j') A (a; + a;') that appear in Eq. (ID32[) remain with the effective cumulants, 
because they appear in the contributions of all the proper diagrams of G^ii'^^ii ^^)- then have 

A{u'-u,) 5(j.Ji)<5(j.J')A(c.'-L.i) Ml[Ui,,-w^-u,-i',uj\u') 

XA (u'l - U2) 5 (j^,j2)'5 (j,;,j'i) A {uj[ - W2) Ml{l,^{h, "W2, -M2;ji,^i,Mi) 

X A(4-U3) (5(j,;j3)(5(j,J^)A(a;2-a;3) Af^^^^ (jg, -wg, -7/3; j'2, ^3, M2) 

X 5 {ii:in) 5 (^^,^'n-l) ^{^'n-l - ^n) ^ (j„ , -Wn , -U„ ; j^,. 1 , W^^.l , <_ 1 ) 

X A« + u) <5(j,J)5(j,j:jA(c.; +c.) Mli[{lu:,u-t,^'n,<). (D33) 

(b) The contribution of the n cumulants of conduction electrons 

■ (5(ki,ki)J(ui,Ui)5(CTi,cr'i)(5(wi, w'l) 

lUJl + Ui £ (ki, CTl) 



X ■ (5(k„,kJ,)(5(u„,0(5(cr„,cr^j)J(tJ„,a;^,). (D34) 

lUJn +Un e(k„,(T„) 

(c) The contribution of the 2n interaction edges 

«(ji,ai,ki,cri,Mi)w(ji,ai,ki,cri,-'Ui) 

X On' "n' o-^j, u^) V (j„, a„, k„, cr„, -Un) ■ (D35) 

As in Eq. (|D4|) there is a factor (—1) for each interaction parametei'^ v (js, as, kg, cts, Us),and they all cancel 
out in pairs for the C/^^, (k, icj), and the C/^^, (k, iw^), but for C/^-^, (k, iws), and tj^^, (k, icjs)^ne of these 
factors (—1) remains and a change of sign is necessary. This sign is not necessary in Eq. (ID33P because is 
cancels out like in the QI/^, (k, iuj). 

(d) The two external lines correspond to the same site ji because this is the only site with local / states, so that 
we have (5(j, j J^(j', j,). 

(e) A factor ±1 obtained employing the rules given in Appendix|Xl and the same arguments used in Section [D II 
can be applied here. In particular, we assume that the contributions to the effective cumulants M'^^J, {],ujs) 
have been already calculated with their correct sign. The factor is therefore 



(D36) 
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FIG. 3: GF diagrams in real space and imaginary frequency of the PAM, with n+1 effective cumulants: we write j', a' , u' , oj' = 
jo, QfQ, Uq, u)'(, and j, a,u,u) = j„+i, a^+i, u„+i, tDn+i. The same diagrams describe a single Anderson impurity (SIAM) at site ji 
when there are local / states only at this site. 

(f) A factor 1/g calculated from appendix IB] We assume that all the factors that appear in the contribution 
of all the diagrams corresponding to the effective cumulants have been already included in those cumulants, 
and we then only need to calculate the g corresponding to a chain, that is g = 1. We can then write 

- = 1. (D37) 
9 

3. Sum the resulting product with respect to 

(a) The site labels js of all the FV for s = 0, 1, . . . , n, but all these sums disappear because there is only one site ji 
with local / states: 



<5(jiJs+i)'5(j«Js) 



(D38) 



with = j' and j„+i = j. 

(b) The momenta ks, the frequencies Ws and the indices Ug of all the internal edges. Because of the delta functions 
in Eqs. (|D33|D34| we have 



U — Us — ^ —u. 



(D39) 



and 



(D40) 



(D41) 



for s = 0, 1, 2, . . . , n. 

Notice that for real space, the sum over momenta does not reduce to a single term, but there is a summation left 
at each CV, and we shall discuss this summation at a later stage, because we have to consider the dependence 
with ks of the factors in Eqs. (|D34ID35|) . 

(c) We have the sum over all as, a'j,, for s = 1, 2, . . . , n and we shall use matrix notation to simplify the calculation.. 

(d) Because of the 5{(Js,(j'g) in Eq. (jD34[) and the spin conservation in the effective cumulants when [az]T-C\ — 
(even if a'^ ^ cks+i) there are no sums left over the labels ct^. 

We shall now consider the contribution of the factors in Eq. (jD33|) : employing Eqs. (jD39IID4T|) we can write 



(D42) 



-ui;y,u}', u') = Ml[l,{]i,uj, u; ji, -w, -m). 



(D43) 
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and 

Maii,a'^Os+i, -i^s+i, -Us+i;3s,uj's,u',) = M^f;^^^, (D44) 

for s = 1, 2, . . . , n — 1. If we now define 

M^^f(j,z^,u) = ArJJij,uj,u;j,^Lu,~u) (D45) 

(recall that in Eq. (|D32p it is u = —1) we can write the contribution of the n + 1 factors in Eq. (jD33l) in the following 
form : 

All these factors are independent of all the and k^, and can therefore be factored out from the summations over 
these wave vectors. 

We still have to include all the factors from Eqs. (|D34ID35P in the contribution of the diagram with n + 1 effective 
cumulants. We calculate v (js, a^, ks, Us, —Us) v (j'^, a's, k'^, cr^, u^) employing Eqs. ()27ID39|D40p and the conservation 

of (Is = 

v{ys,ct's,K,<^s,u'Jv{js,a,,'ks,(Js,-u,) = ■^V{a'^,ks,<T,u'^)V{as,ks,a,-Us). (D47) 
Combining Eqs. (|D34|D35|D47p with Eqs. (|D15IID20|D23|) we write 



Wa'.a (k, a, u' = +; iws) = i^^* {a', k, a)V{a, k, a) 



Ns ' ' ' ' + e (k, a) 

= ^V* (a', k, a)V{a, k, (k, -lUs) (D48) 

in place of Eq. (lD20t : 

We return to the GF Gaai{i, it^s) defined in Eq. (|D32p . which corresponds to u' = +, so that the factors from Eqs. 
()D34|D35P could be put then in the form (cf. the labels in figure [3]) 

cr, +; iujn-i) ■ ■ ■ Wa' a2 (k2, o", +',1^2) Wa' ai (ki, cr, , (D49) 

and still have to be summed over all the ks. Employing Eqs. (|D23IID26|) we introduce 

* k 

where we have taken the sum over k, used again Us — — w, and employed z in place of the Matsubara frequency iuj. 
The sum over all ks of the contribution in Eq. (jD49|) then becomes 

Wa'^,a„ (cr; z) Wa'^_^,a^_, (ct; z) . . . Wa'^,a2 (c^; z) Wa'^^a^ (o"; z) . (D51) 

To simplify the calculation we now introduce again two matrices (cf. Section [IV 1[) 

{M}^^^^ArJJ{j,,z,u), (D52) 

and 

{W}^^ ^^Wa'^a{'J,z). (D53) 

The contribution of the diagram with n+1 cumulants takes then the same form of Eq. (|D30p , but with Eqs. (|D52ID53p 
in place of Eqs. (|D28ID29p : 

{M}_, {W}^, {M}„„„. _^ {W},, . . . {W}„, {M}„^„, {M}„^„, = 

{(M . W)" . M}„_„, = {M . (W . M)"}„_„, (D54) 
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Appendix E: Calculation of approximate GF 

Here we define and give the formal expression in the Atomic Approximation of the GF for the PAM and the SI AM, 
that were left out in Section IVl CI 

1. The other approximate GF for the PAM 

a. In reciprocal space and imaginary frequencies 
The approximate GF Qj^^, (k, iuj) in reciprocal space and imaginary time is defined by 

((r(/;k,a,w = -,c^) Y{c-y,a',u',u'))^)=gil,{\<i,iu:)^{u + u') ^{uj + uj') 5(k,k'), (El) 
and employing the Rule 3.7 in Section Hi Al we obtain 

gil, (k, iu^)^-Y. Gil (k, tuj) V{a\ k,a',u = +) (k, zco) , (E2) 

a' 

where 

V{a',k,a',u = +) ^V*ia',k,a'). (E3) 
We now introduce a column vector G^f (k, iuj) so that 



where we have changed the dummy variable cr' into a, and we must remember that Gq^ g = G^^ ^ = 0- 
Substituting Eqs. (|E3l72|73p into Eq. we obtain 

ro n ■ \ I ™ii + ™13 1 

G('=''^^(k, ico) = -V* ^""'^ '"''^ (E5) 

1 - \V\ G^^^ (k, luj) {mil + + mis + m^i) 



Gl''^P{k, iw) = -V* 2 ^ ^ (E6) 

1 - \V\ Gc^i (k, luj) (m22 + - TO24 - TO42) 

We define the approximate Ga{j' O^j *'^) with 

((r(c;k,a,w = -,c.) Y{f-y,a',u\io'))^)=G:i,{k,iiu) A{u + u') A(c. + c.') (5(k,k'), (E7) 
and from the Rule 3.7 in Section Hi Al we obtain 

G^l (k, iio) = -J2 5t (k, V{ai ,k,a,u = -)Gil^, (k, zu), (E8) 

ai 

where 

V{ai,k, a,u = -) = V{ai,k, a). (E9) 
We now introduce a row vector G'^J''^^lk,iuj) so that {GJ^'^'''P(k, io;)}^, — (k, iw), and then 

G^/'^^'(k, ^u^) ^ (GlUk, i^) , tJ:^_.,(k, i^)) . (ElO) 
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Substituting Eqs. (|E9I72I73I) into Eq. dES]) we obtain 



Gf''P{k,iLu) = -V __ IV ^g^j) 

1 - \V\ g°,. (k, iuj) (mil + 77733 + wi3 + 77731) 



, GciO^^"^^) ('^22 - ,77724 - 77744 

Gf '"f (k, iuj) = ± ^ (E12) 

1 - \V\ CJ"^ (k, iu) (77722 + "144 - '7^24 " nT-42) 

Finally, we define the approximate Q^'^(k,iLu) with 

((y(c;k,(7,t7 = -,w) y(c;k',a',7i',cj'))+) =e^"(k,iw) A(77 + 7i') A(cj + cj') (5(k',k) 5(a,a'), (E13) 
and using Rule 3.7 in Section Hi Ap we obtain 

1+ n«l,k,a,77 = -)e;|;f„,(k,^c.)T/(a;,k^a,77 = +)^,%(k',^^) U(k,k'). (E14) 

We also introduce the scalar G^'^'"P(k, iw) so that 

G^^^'^f (k, iij) = ^^(k, k, iuj). (E15) 
Substituting Eqs. (|E3|E9l72l73p into Eq. (|ET41) we obtain 

t^ccapf, . . ro fi. ■ ^ , \V\^gl^{Ki^) {mil + 77133 + mi3+m3i) . 

G^ ''(k,7a;) = g^^^{k,iuj) + - . n ■ ^/ T^^ct (k, *'^) 

1-|V| t/^' .J. (k,7i:jj (7nii + 77733 + 77713 + 77731) 



Gl^ (k,iw) 



1 - I g'^.^ (k, 7i:j) (t^h + 7^33 + 77113 + 77131) ' 



and 



I V^l^ g°^ (k, iuj) (77722 + "144 - "124 - "142) ^0 



G ^'(k,7w) - g,^^(k,7a.) + - r^c,;(k,»w) 

1 - I I g^^l (k, 7W ) (77722 + m-44 - "l24 - "^42 ) 



1 - \V\^ (k, iuj) (77722 + "144 - "124 - "142) ' 



h. In real space and imagtnary 



(E16) 



(E17) 



We follow the same procedure used in section IVIB 21 to derive the GF in real space when the / electron is created 
and destroyed at the same site. Considering again a rectangular conduction band we find G^^(7aj) by integrating Eqs. 
(111161): 

G/c,ap(,^) = ( ^ot.T(*^)^ ^ V ('^''^) ( mil + mi3\ 

^ \gil^{lUj) ) V 1 - 1^1' ^c.T ^31 +"^33^' 



so that 



G^'^^M = In f 4#4^^) r " + , (E18) 



21? - - M/ \ "131 + 77733 / ' 
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where 



and in the same way we obtain 



V* f AcriiUj) + D - fl\ / 77122 " "^24 

In ' ' ' 



In a similar way we obtain the G^^^°p( 



2D \Af^{iuj) - D - fi / \ m42 - m44 
(oi'^oAi^) ^Si%Aioj)) by integrating Eqs. (|Ell | E12p : 



■ In ' 



2D \A^{iuj) - D - 



(mil 



TO31 , mi3 + TO33 



and 



pc/,ap,. s V f A^jlLj) + D - fj. 



77122 - "^42 , "724 - "l44 



To obtain G^'=(iw) we integrate Eqs. (|E16IE17|) : 

Employing this relation, we can then write 

G("'"P(iw) = -V* G^ 



1 / A^jiuj) + D - fi 
2D ^ \Aaiiio) -D- fi 



llUJ 



mil + mi3 
77731 + 7M33 



(E19) 



(E20) 



(E21) 



(E22) 



(E23) 



(E24) 



Gf '"PCiw) = -y* Gf '"^(iw) ' ""22 ^24 



77742 ^ '7744 



G^^'"P(ia;) = -y G^'^'"^(ia;) (77711+77731 ,77713 + 77133], 



Gl^'°'^{iuj) = -V Gf '"^(iw) ( 77122 - , ^24 - TO44 



Gf = M7 + |y|'Gf(iw) 



Gf M^Mf + l^rGfM 



1 1 



1 1, 



(E25) 

(E26) 
(E27) 

(E28) 
(E29) 



c. Green's functions with the usual Fermi operators f and . 

It is interesting to calculate the Gf of the usual Fermi operators, related to the Hubbard operators through 

f^Xoa + (T X.,d (E30) 
where for typographical convenienc we use —a in place of a. It is straightforward to obtain 



((/.;/t))^(k) = 



(E31) 
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where we used Eqs. (|94l95p . 
In a symilar way we find 



((c,(k);/t)) =-V QcA^^'^)^l^ (E33) 



and 



((c.(k);4(k'))) =<5(k,k') 5" ^^''^^ jj. (E34) 

As in the previous section we calculate the GF in imaginary frequency and real space when the / electron is created 
and destroyed at the same site 

((/.; /t))^ = MfJ [l + \V\^ MffGl^{^Lo) 
We also obtain 

{{f.;cl))^^-V* Gl^{rio)MfJ, (E35) 

and 

{{c.■Jl))^^~VGf{^^)Mff. (E36) 
The ((c.;4)), is given by Gf{z) in Eq. (|E23| 

2. The other approximate GF for the SIAM 

a. With f electrons in real space, and imaginary frequencies 

The approximate GF Q^'^,{ji,'k' ,iuj) with the impurity at ji is defined by 

{{Y{f;i,a,u = -,u) Y {c;k' ,a' ,u' ,u:'))^) ^ G^.i^^ ,icj) A {u + u') A + w') <5(jJ,), (E38) 
and employing the Rule 3.7a in Section Hi Bl we obtain 

Gil, a« = 0, k, iu:) = -J2 S"' (j^ = 0' "(j =0,a\k,a\u = +) Q%. (k, lu:) , (E39) 

a' 

where 

v{i = 0,a',k,cr',u = +) = N7^V*ia',k,a'). (E40) 
We now introduce a column vector G;Jf (j; = 0, k, iuj) so that 



(E41) 



where we have changed the dummy variable <t' into a, and we must remember that Qq^ ^ — qI"^ g — ^ 
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Substituting Eqs. (IE40I72I73|) into Eq. (|E39)) we obtain 

mil + "^l3 



i^Gc -f (k, iujj 2 

'Ns ' 1 - 1^1 <y5t(*'^) C*^!! + "^33 + "^13 + JTlsi) 



G(='"^'(k, 7LU) = -^Gl^ (k, ZW) 2 . ^- . (E42) 



/ 77122 - '7124 



Gf ■'^^'(k, iij) = -^g°^ (k, icv) 2 ^ ^ (E43) 

vNs 1 - 1^1 (»7l22 + W44 - TO24 - W42) 



We define the approximate tj^'^, (k, j^, itj) with 

{{Y{c;k,<j,u = -,u) Y{f;y,a',u',io'))^)^g^i,{Ky,7u) A{u + u') A{lo + lo') 5(j'JJ, (E44) 
and from the Rule 3.7a in Section Hi Bl we obtain 

g^l (k, j,; =0,iu:) = -Y, Gla (k, v{i = 0, ai, k, a, u = ~)Qi{^^, (j,; - 0, iu:), (E45) 

where 

=0,ai,k,a,u = -) =iV7V(ai,k,cr). (E46) 
We now introduce a row vector GJj-''''^^(k, j,; = 0, iw) so that {GJ^'^'°P(k, = 0,«cj)}^, = ^^{./(kj^ = 0,icj), and then 

G^/''^f(k,j. = 0,zw) = (t;:^o,(kJ, = 0,1c.) ,C/:^_,,(kJ,; = QM) ■ (E47) 
Substituting Eqs. (|E40|72|73p into Eq. (fEiS]) we obtain 



V (mil +"731 ,™i3 + m33) 



Gf '"P(k, iw) = (k, iLo) 2^^— " ' "-^ (E48) 

V^s ' 1 - 1^1 {mil + 7n33 + 7ni3 + Tn^i) 

y (to22-TO42 ,77124 -"744) 

Gf -(k,!..) = -^Ql, (k,7-) J 1 (E49) 

V^^s l"!*^! </'i(7W) (77722 + 77144 — m24 — 77142) 

Finally, we define the approximate t/^'^(k, k', ilj) with 

((y(c;k,a,u = -,c.) y(c;k',a',u',L.'))+) =e^'(k,k',ic.) A(u + 7i') A + w') J(j',j,) 5 (a, a'), (E50) 
and using Rule 3.7a in Section HlBj) we obtain 
gr(k,k',ic.) = g,% (k,ic.)x 

,5(k,k')+ 7;(j=0,ai,k,a,u = -)g^f^,(j,ic.)i;(j = 0,a'i,k',a,u-+)g,%(k',H [• (E51) 

We also introduce the scalar G^'^'"P(k, k', jcj) so that 

G:^^''^P(k, k', iu) = 5^=(k, k', iu). (E52) 
Substituting Eqs. (|E40l72l73p into Eq. ((E5T|) we obtain 

.^0 /, ■ N (77711 + "733 + 77113 + 77I31) ^0 /I / • \ ^T^roN 

G^ ''(k,k ,iw) = e;^.^(k,iw)5(k,k) + — — g;^.^(k,ia;) —2 — — -gi^{k,iuj) (E53) 

1 - (^tl*"^) + "^33 + 77113 + mai) 

Gr-^(k, k', 7.) = Gl, (k, ^.) 5 (k, k') (k, 1.) ("^-22 + 77744 -7.24 -"742) ^^o^ ^^^^^ 

^^s <i5;(7'^) (77722 + "744 - 77124 - 77142) 
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b. Green's functions with the conduction electron in the Wannier representation. 

In the impurity case, it is more convenient to use the GF with the conduction electrons in the Wannier representation, 
locahzed at the impurity site. To that purpose, we employ Eq. 



^1- = 4f E '^^P(-*k . R,)Cl , ((->23)) 



and as before we use Eqs. (|94l95p : mI-^ = (tou + TO33 + mi3 + TO31), and Alf-^ — (to22 + W44 — ^7124 — ^42) . 



— V'""!! < '"-33 I '"'in ! "I'iLJj cliiv-i i>-i^ 

We now apply Eq. ([231) to Eq. ((E42)) with R,; = 0, and find 



mil + 



Gf = "^o;/ = ^ct E55 



Employing Eq. ((79|) we obtain 



Gf-^(za;) = -y* ^ [mil +mi3\ ^^56) 

^ 1 - |l^r ni^uj)M(^ \m3i + m33 y 



and in a symilar way we find 

GfrP{^u) ^ -V* ^ "^A . (E57) 

^ 1 - \Vf ipi{iuj)Ml^ \m42 - m44 / 

To obtain the G^-'^'°P(iw) = , ^/f employ Cja = X^k *5^P(+*'* ' R-jOC'kfT, and in a similar 

way we find for = 

;f --(zc.) = -V ^ 



G^-^'"P(ia;) = -F;; (tou +m3i ,77113 + 777.33), (E58) 

and 

Gf'^Piiw) = -V- 1^^!^'^^'')'.^ ,,^,// ( "122 - m42 , TO24 - 77744 ) • (E59) 

1 - 1^1 '^;(«w)Mj;-' 

The remaining G^^'°^'(7cj) are now easily obtained 



G-'«''(7c^) = jj. (E60) 



Employing this equation we can now write 



G 



Gf •'^P(7w) = -y* Gl'^^'Piiuj) i + '^13 ) , (E61) 

I \TO31+m33/ 



Gf •"^(iw) = -V* Gf '"^'(iw) ""22 m24 ^ ^gg2^ 

^ ^ \r7742 - TO44 J 

'"^(iuj) = -V G^'=^°^(ia;) (77711 + 7^31 , 77713 + "133) > (E63) 



Gl^'°'^{iuj) = -V G!^'''''^(7a;) (^77122 - 77742 , "724 - 77744) 1 (E64) 
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and also 



/ mil mi3 \ 

^t^'"'(*^) = . i.r?' "r .rff " 1^1' ("^""^=^3 - miamai) ( \ ' A , (E65) 

1 — |l/| iy9^(iaj)Ai| ^ y— i i y 



TO22 ^24 
17142 m44 t 

'^V'^'i'^) - , ,2 // - Gr(*^) ('^22m44 - m24m42) I ; ; 1 . (E66) 

1 - 1^1 ^i[ii^)M^-' \l 1 

We can now rewrite these two equation in the form (cf. Eq. ([83 



Gf/'"P(ic.) = + \V\^ Gf (*c.) I M^f - (mnm33 - migmai) f M , (E67) 



Gf^''^P(ic^) = + 1^1" Gf (zc^) |Af/-^ - (m22m44 - m24m42) j j | ■ (E68) 

c. Green's functions with the usual Fermi operators f and . 

As with the PAM, we calculate the Gf of the usual Fermi operators, (cf. Eq. (IE30|) ). It is straightforward to obtain 

M-ff 

1 - \V\ ipa[z) Ma 

where we used Eqs. (|94l95p . In a symilar way we find 

V* Mff 
((/.;ct (k)))^ = --=t;°.(k,z) 2 \ , (E70) 

V 1 - |y| (y5cr(z) il/^-' 

and 

((c.(k); ct (k'))), - (k, zc.) S (k, k') + (k, z) - (k', z) . (E72) 

^Vs ■ 1 - \V\ (Pa{z) My 

As in the previous section we calculate the GF with the conduction electron in the Wannier representation 

[fy^4=o.)) = ^V* Gl^iiLo) Mff , (E73) 



and 

((cj.o.; /J)), = 'V Gl^izuj) Mff. (E74) 
The ((cj=oa;c]^o.))^ is given by G^^(iw) in Eq. (jE60| 

l-\V\'^,iz)Ml 
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3. Summary of the approximate GF for the PAM 

a. GF in reciprocal space and imaginary frequency 



GV'-iK i^) = , „ ' 'I , (E76) 



and 

/ m22 

.ff.ap- - ■ V"^42 m44 



\^\^ Qc,i (k, Z) (m22TO44 - m24TO42) ^| j 

1 - \V\' (k, 2;) (m22 + m44 - m24 - m42) 



For 

7/^ 



(k, iw) 



'mil +mi3 \ 

I 77131 + m33 j 

-V 

1 - \vf' Q^.^ (k, iuj) (mil + ^33 + mi3 + m3i) 



(E78) 



Gf'"^(k,ia;) 



-V* 



m22 - TO24 

, m42 - m44 , 



1 _ gO^ (-Jj^ (-^22 + m44 - m24 - m42) 

For G^/'«f(k,ia;) = (gi%A^M ,<'=^rf(k,iu;)): 



Gf'"f(k,iu;) = -V 

1 - \V\ (k, iuj) (mil + m33 + mis + »Ti3i) 

QcS (k, iuj) [17122 - m42 , m2i - m44 ) 

Gf '"P(k,iu;) = -V ^— ^ ^ 

1 - g°^i (k, iw) (m22 + m44 - m24 - m42) 



(E79) 



(E80) 



(E81) 



1 - \ vf Q^.^ (k, iw) (mil + mss + mi3 + msi) 



(E82) 



Gr'"^'(k,zw) = 



Oil (k,ic.) 



1 - 5°4^ (k, ico) (m22 + m44 - m24 - m42) 



(E83) 
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64. = m22m33 - m24m42. 



39 



c. Green's functions with the usual Fermi operators f and p 



«/.;/t)>^(k) = 



(E91) 



((/.;4(k))) =-y* 



(E92) 



((c.(k);/t)) =-y 



(k,z)M// 



(E93) 



((c<.(k);ct(k'))) =Mk,k') 



l-\V\' GlAKz) Mlf 
In imaginary frequency and real space when the / electron is created and destroyed at the same site 

((A;4(k)))^ = -y* G-MM//. 
{{cy,fl))^ = -V G-MM//. 



(E94) 

(E95) 

(E96) 
(E97) 



^ 2D \A„{z)-D-n 



(E98) 



4. Summary of the approximate GF for the SIAM 

a. GF with conduction electrons in k space 



Qff,ap 



mil mi3 \ ,,,,2 /• X / ^ I 1 ^1 

\-\V\ (miim33 - mismsi) 

^TO31 mSS/ 1 ^ 

1 - ¥'t(^'^) + '^33 + mi3 + msi) 



(E99) 



Qff,a.p 

4- 



- |F| (^4.(1^) (m22m44 - m24m42) ^ ^ 

^TO42 m44 J \^ i- j 

1 - \Vf' ipiiiuj) (m22 + m44 - m24 - ^42) 



(ElOO) 



Gf (k,iw) 



V* 

7W. 



mil + mis 
mai +m33^ 



1 - |y| (/?t(*'^) + '^ss + mis + msi) 



(ElOl) 
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/ TO22 - TO24 \ 

. T/* \ 77142 — rn44 I 

Gf (k, ia;) = -^Ql, (k, ia;) 2 — ^ : (E102) 

, V n (mil +77231 ,TOi3 + 77133) 

Gf(k,iu;) = =e?o^(k,iu;) 5^ ^ (E103) 

V-'Vs ' <^t(*'^) ('^ll +^^33 +mi3 + m3i) 



„ F „ I ™22 - m42 , 77124 - TO44 ) 

Gf (k, i^) = -^Gl, (k, i^) } 1 (E104) 

Gf (k, k', ^.) = Gl, (k, ^u:) 6 (k, k') + ^g", (k, ^u:) (^11 + m33 + mi3 + m3i) ^ (j,', (E105) 

-'Vs ¥'t(*'^)(™ii +"^33 + ^13 + 77131) 

Gf (k, k', ^.) = Gl, (k, ^o:) 5 (k, k') + 1^^°, (k, zu,) ^ i-nn,, ^ n:,^ n.,, - rn,,) 

^ ^ ' Ns 1 -</94.(^a;) (77122 + 7^44 -m24- 77742) 

6. Green's functions with conduction electron in the Wannier representation. 



G(-^(za;) = -y* jj + ) , (E107) 

^ l-IVlVtM^t V"^3i+rn33y 



Gf ^=(70;) = -V* ^ ~ I • (E108) 

1-\V\^ (pi{ioj)Mlf \m42-m44J 

Gf(7a;) = -y— — ,mi3 + m33), (E109) 

1 — |K I iy9-f (zw)M| ^ ^ ' 
Gf{iuj) = -V- ,^^,2^^'^.^ (m22 - 77242 , 77724 - m44) • (EllO) 

1 - |V"| (^4.(7u;)Mj;'' ^ ^ 

G^^(zw) = ^^^^^ rj. (Elll) 

c. Green's functions with the usual Fermi operators f and f^ . 



((/- fl)). = , |,,|2 ^ , , ,,ff (E112) 

1 - |y| iy9a(-2;) Mi-^ 



/<.; 4=0.))^ = -1^* G:'(«'^) M//, (E113) 
((cj=o.; 4)), = -^^ G-M M/A (E114) 



Appendix F: The free /-electron GF 

We first collect a few equations from the previous appendices (cf. Eqs. (|C1IC4[) ) 



(cf. Eqs. Jell) 



^==\ [YilUn - T2)y(72, 0) j ^F{n- T2) , 



5(71, n; 72, 7-2) = exp {I3Q) Tr{exp {-13%) exp {t{H)Y^^ exp {-t^H) exp {T^l-i)Yy^ exp {-t^T-L)} . 
(cf. Eq. (ICT2I) ) 







^;(7i,wi;72,W2) = A(a;i +W2) / dxF (a;) exp [ixwi] , 

and (cf. Eq. EH 



(f(/; k, a, u = -, 0.)^ (/; k', a', w', ^')) = ^^'(k, A (^i + ^^') A (^.k + u'\^) A (c^ + u') . 

These equations correspond to the exact GF, and for the free GF (i.e. with no hybridization) we shall now 
that 

where 

Substituting the general ^(7, r) operators by X operators we have 



<y M = ^(f(7,r)y(7',0))^^^ ^= ^(x,,„(r)4„,(0)) 

= exp m) Tt {exp (-/3Ho) exp (tHo)X„ (b^,) exp (-rHo)Xj (,,^,) } . 
Now we use a basis {|c)} of eigenstates of "Hq to calculate the trace: 

(^) = exp (/?!]) ^ (c| exp (-/3Ho) exp (tHo)X,- (b,,) exp (-rHo)X,,(,,^b,) |c) 

C 

= exp {pn) ^ (^b'.c^a'^a^cb exp {{t - P) Eb) exp (-tEq') = 

C 

= <5a/^Q(56',b exp (^rj) exp {~l3eb) exp (r (e^ - Sa) ) =(5qq'F^ (t) . 



We calculate the Fourier transform 

r/5 



Gia°{K = (5qq' / (x) exp [iTw] = (5q,q/ exp (/3fi) exp (-^flcf,) / dx exp (a; (e^ - Sa) ) exp [ixw] 

Jo Jo 



6aa' exp (/3rj) exp {-Pet)-. — {exp (/3 {iuj + £„))- 1} , 
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and employing exp (i (3uj) = — 1 we find® 

nff.o,, ^ . exp (/3f^) [exp (-/3£h)+ exp . £>a 

y'^-l^, (k, = - daa' : — == - Oaa 



Da 



+ So 



(FIO) 



where 



lUJ + (Eb - Sa) 

Da = {Xaa + X^h) = exp {(5^) [exp (-/3eb)+ exp (-/3£a)] • (Fll) 



Appendix G: The approximate GF for oo and a rectangular band in the impurity case. 

We assume that the dispersion relation of the conduction electrons corresponds to a rectangular band, with 

-D< Ek., < D. 



(Gl) 



In the present case we have only a single a = (0, cr), and we assume that y(a', k, cr)!^* (a, k, ct) = \V\^, so that from 
Eqs. P23ID50P 



\v\' 



z — e (k, z) 2D J z — X 



dx. 



(G2) 



To avoid singularities we employ z = uj + is with s > 0, so that 



W{uj + is) 



2D 



In 



UJ + is — (D — ii) 



UJ + is — {—D — /Li) 



and (cf. 



A{uj + is) 



\V\ 



In 



UJ -\- is — {D — ^) 



M^f^iuj + is) 



2D [uj + is- {-D- fj.) 
For the band with zero width we take E'k.o- — (cf. Section |Vl, so that from Eq. (IG2 



W'iz) 



-\VY 



z- E^ + fi 



(G3) 



(G4) 



(G5) 



* Rewriting Eq. I)F4|| for the unperturbed case, but with u' — > —u' and ui' — > —uj' (after exchanging primed and unprimed variables), we 
have 

( (y{/; k', a', -u' = -, -J) Yif; ; k, a, u,u;))^)^ = f',"(k', -u;') A (« - u') A («k - «'k') A {uj - u') , 
and employing Eq. I)F10| | we find: 



Da 



Da 



-iul' + (Cb - Sa) il^' + {Sa - Efc) 

For the conduction electrons we should then have (putting Da = 1 and (sa — £{,) = e (ki, cri) — = e (ki, (Ti)) 

((C(k'i,<7i,-ui = C(ki,<7i,«i,a;i))^)^ = gS'='0(k'i,-u;i) A («i - <) A («iki - ^ki) A - 

5(ki , k'l )(5 (mi , m'i ) (5((Ti , o-i ) 5(tJi , a;i ) , 



+ e (ki, (Ti) 



which is just Eqs. I|30I37| I for u\ = +. Note that one then has 
and (cf. Eq. (g2j) 



-1 



4-e(k;,a0 



iw' - (Sa - St) 



which are used in another context. 



43 



We then have from Eq. (lG4t (cf. also Eq. (1591) ): 



M"'' (z) = (l + G^^^°* (z) • W"*) • G^f-"' (z) = 



{z- ES + /i) G-^'^^'^* (z) 



and substituting Eq. (|G6]) in Eq. (lG4l) (cf. Eq. ([55])) we find 



A° 



2D 



In 



a; + is — (Z? — /i) 



(z - + /i) G^-'^-''* (z) 



+ is-{-D-^I)\ z - E'^^ + ^1 ~ \V\^ G//>* (z) ■ 



Employing Eq. (|57p we can now write the approximate GF as 
G^f = M (I- A)"^ = 



(z - E- + m; 


G//.'^* (z) 


(z - + /i) - 11^1" G//.-* (z) + ^ In 


u+is — (i) — 


(z - ES + fi) Gff'^t (z) 


LJ + 'is — ( — Z) — /J-) 



(G6) 



(G7) 



(G8) 



Employing the GF of the PAM with a band of zeroth width to calculate the impurity M"'' (z) . 



For the PAM we have 



M'^P (z) =[1 + G^f'"* (z) ■ W j • G-^^^°* (z) . (G9) 

where 

{W}„,„ = W^,^ (k, a, z) = V^, (k, a) V* (k, cj) G° , (k, z) . (GIO) 
Employing the conduction electron free GF for a zeroth width band 

(k,z) = ^4 T = — , (Gil) 

"'"^ ' ' z-£(k,z) z-E§+n' ^ ' 

and assuming that Va' (k, cr) V* (k, ct) = we obtain 



- \V\ 

(G12) 



which is exactly what we obtained in Eq. (jGSp . 
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14) = \d,\.) 
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15) = |d,t) 
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16) = |d,n) 



TABLE I: The table gives the matrix elements of the six operators Xo+ = X+ = Yi, Xo- = X- = Vi, X^ii = T+ = Yij, Xo-d = 
r_ = Ki, = Ys, and = Ye in the basis of the sixteen states defined in the last column. The matrix is separated into the 
sub matrices (n| Y, \n') connecting states with n = 0, 1, 2, 3 electrons to states with n' = n + 1 = 1, 2, 3, 4. The value of the 
matrix elements is either 1 or —1, as indicated in the table. We use \d) = \-\ — ), as in Table HIl to emphasize that X± \d) = 0. 
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|2) 










-|5) 


5 




10) = |+,t) 


|3) 
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11) = \d,0} 
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12) = \-,n) 




|6) 






-|7) 


|8) 
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13) = i+,t;) 
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|10) 
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111) 
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15) = |d,t) 
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|10) 


111) 




1 




16) = \d,U) 






|12) 


|13) 


|14) 


-|15) 



TABLE II: The elements in the table give the state that is the result of applying the destruction operators on the top of the 
table to each of the sixteen states defined in the second column, where we use |d) — | + — ) to indicate the state with two local 
electrons (note that Xo+ \d) = Xo+ | + — ) ~ and it is neither |— ) nor |0)). The numbers in the first column gives the ordering 
of the states in the local subspaces with n = 0,1, 2, 3, 4 electrons 



